LINEAR  SI6NAL  PROCESSING  AND  ULTRASONIC  TRANSVERSAL  FILTERS. (U) 


AD.  AO  6?  0 45 


■'.ST?,-  jfcV/^  *■»**:' -•  ?>,;^»-v;.‘ci  ’ 


J 


A. 

§ 


L 


1 m 


ADVANCE  COPY 
LEVEL 

" '*r*>V  *9  Bfc,  _ 

> r,  n*- 
* i-  vn  . ..  u pj?  & 


l&'Mosr  P JJecf  -3 

(P 


NAVAL  UNDERSEA  RESEARCH  AND  DEVELOPMENT  CENTER 


Pasadena,  California 


(JU. 


Nov< 


Microwave  Theory  and  Techniques , 
November  1969. 


DISTRIBUTION  STATEMENT~A~ 


Approved  for  public  release; 
Distribution  Unlimited 


*1 


ABSTRACT 

The  role  of  linear  transversal  filters  in  signal  processing  is  discussed  in 
Section  I.  Linear  filters  for  signal  processing  must  often  have  complicated  im- 
pulse responses,  with  large  bandwidth  and  large  time-bandwidth  product.  The 
linear  transversal  filter,  a delay  line  with  weighted  and  summed  taps,  is  ideally 
suited  for  the  implementation  of  such  filters  because  of  its  simplicity  of  synthe- 
sis. The  filter's  impulse  response  is  derived  by  the  application  of  some  con- 
cepts from  the  theories  of  vector  spaces  and  sampling,  and  is  shown  to  be  equal 
to  the  tap  weighting  function.  Thus,  the  synthesis  procedure  consists  merely  of 
sampling  the  specified  impulse  response  at  appropriate-intervals  and  using  the 

sample  values  as  the  tap  weights. y m.t< 

The  utility  of  the  transversal  filter  in  signal  processing  is  illustrated  by  an 
example  from  scatterer  distribution  mapping.  The  illustration  is  applied  to  two 
hypothetical  systems  — a sonar  and  an  astronomical  radar.  In  both  these  cases, 
it  is  not  possible  for  a single  filter  to  process  the  signal  in  real  time.  Signal 
processing  in  compressed  time  is  discussed  as  an  alternative  to  the  use  of  a 
large  number  of  filters  in  parallel.  If  the  processing  filter  has  a bandwidth  capa- 
bility in  excess  of  the  signal's  bandwidth,  the  signal  can  be  time  compressed  and 
processed  serially  in  time.  A generalized  receiver,  employing  time  compression, 
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frequency  translation  and  multiple-output-port  transversal  filtering  is  developed 
from  these  ideas. 

In  Section  n,  a generalized  transversal  filter  is  described  and  analyzed.  A 
delay  line  with  multiple  arrays  of  taps,  each  array  with  a multiplicity  of  weight- 

I I 

ing  functions,  has  as  the  impulse  response  between  any  pair  of  ports  the  cross- 
correlation function  of  the  weighting  functions  for  the  two  ports. 

A number  of  implementations  of  transversal  filters  employing  a variety  of 


delay  line  types  are  described  and  some  aspects  of  transduction  and  wave  propa- 


INTRODUCTION 


Linear  systems  play  a leading  role  in  science  and  engineering.  In  particular, 
linear  filters,  which  occupy  a position  of  universal  importance  in  the  electrical, 
electi’onic,  and  communication  sciences  have  been  studied  extensively.  There 
are  three  major  classes  of  linear  filters.  The  earliest  developed  was  the 
"longitudinal"  filter,  in  which  the  signal  is  operated  upon  successively  by  the 
structural  elements  of  the  filter.  The  ladder  netwoi’k  and  the  meander  line  are 
typical  examples.  Because  of  its  importance,  and  since  it  is  the  earliest  of  the 
three  types,  the  longitudinal  filter  is  in  an  advanced  state  of  development.  The 
second  class  is  the  "recursive"  filter,  in  which  previous  values  of  the  signal  are 
modified  and  recombined  with  new  values  through  feedback  loops.  The  linear 
feedback  network  is  a typical  example  of  the  recursive  filter,  which  also  has 
been  developed  to  an  advanced  state. 

The  third  and  newest  class  is  the  "transversal"  filter,  in  which  the  signal 
propagates  in  a delay  line  with  transverse  taps  distributed  along  its  length:  the 
output  is  formed  by  weighting  and  summing  the  contributions  from  the  taps 
(Fig.  1).  The  transversal  filter,  originally  conceived  for  the  equalization  of 
television  signals,  has  been  used  extensively  as  an  equalizing  filter  in  telephone 
systems.  In  these  applications,  electromagnetic  delay  lines,  both  distributed 
and  lumped,  have  been  employed,  but  the  total  delay  time  available  in  filters  with 
practical  physical  dimensions  has  been  limited  by  the  large  propagation  velocity 
of  electromagnetic  waves. 


The  use  of  ultrasonic  delay  lines  results  in  a reduction  in  size,  since  the 
propagation  velocity  for  acoustic  waves  in  solids  is  104  to  10s  times  slower  than 
for  electromagnetic  waves.  With  the  variety  of  methods  now  available  for  tapping 
microwave  acoustic  devices,  the  tapped  ultrasonic  delay  line  provides  a practical 
means  for  implementing  transversal  filters.  As  a consequence,  the  transversal 
filter  has  many  interesting  and  useful  applications  in  microwave  technology,  as 
well  as  in  filtering  technology  in  general. 

This  paper  presents  the  theory  and  properties  of  transversal  filters,  along 
with  selected  ultrasonic  delay  line  implementations.  It  will  be  shown  that  the 
transversal  filter,  especially  as  embodied  in  contiguously  tapped  ultrasonic  delay 
lines,  is  an  ideal  means  of  implementing  linear  signal  processors.  The  paper  is 
in  two  sections:  In  the  first  section,  a time-space  equivalence,  based  on  some 
concepts  from  vector  spaces  and  sampling  theory,  is  applied  to  transversal  filter 
analysis  and  synthesis,  and  the  use  of  the  transversal  filter  in  a particular 
application  - the  mapping  of  a distribution  of  scatterers  - is  described  in  detail; 
the  second  section  presents  the  theory  of  ultrasonic  delay-line  transversal  filters 
and  describes  some  implementations  by  means  of  magnetostrictive,  piezoelectric, 
and  optoncoustic  taps. 


LINEAR  SIGNAL  PROCESSING 

Physical  signals  can  be  represented  as  vector  elements  in  the  Hilbert  space 
of  square  intcgrable  functions.  A linear  filter  is  a linear  operator  over  this  space 
and  its  impulse  response  is  the  kernel  of  the  operator.  If  the  filter  is  invariant 
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under  translation,  its  output  is  the  convolution  of  the  input  signal  with  the  impulse 

t * 

response,  which  is  also  an  element  of  the  space.  The  output  also  can  be  repre- 
sented at  each  point  as  the  inner  product  of  the  signal  vector  with  a reversed  and 
shifted  replica  of  the  impulse  response  vector  [1]. 

A signal  can  be  assumed  without  loss  of  generality  to  have  finite  bandwidth 
in  any  physical  system,  and  can  be  represented  by  samples  taken  at  intervals 
approximately  equal  to  the  reciprocal  of  this  bandwidth.  The  sequence  of  these 
samples  is  the  sampled  signal;  the  original  signal  is  reconstructed  from  it  by 
appropriate  linear  filtering  [2],  The  space  of  sampled  signals  is  a linear  vector 
space  whose  vectors  have  the  sample  values  as  components  [1].  By  appropriately 

[■  1 ' 

defining  the  signal's  bandwidth,  W,  and  its  length,  T,  it  is  possible  to  approximate 

the  sampled  signal  vector  by  a finite  number,  TW,  of  components  [3]. 

All  of  the  preceding  remarks  about  sampled  signals  apply  equally  to  signals 

in  time  and  in  space  [4].  The  taps  on  a delay  line  perform  a spatial  sampling  of 

the  signal  in  the  delay  line  and  if  the  taps  are  ideal  (i.e. , if  their  impulse  responses 

are  6-functions),  the  sequence  of  tap  outputs  is  the  spatially  sampled  signal.  If 

the  tap  outputs  are  weighted  and  summed,  the  summation  yields  a linear  operator 

whose  kernel  is  the  weighting  function.  The  sum  can  be  regarded  as  the  output 

of  either  a space-domain  or  a time-domain  filter.  This  kind  of  linear  filter  is 

called  a transversal  filter.  Kallmann  introduced  the  transversal  filter  in  1940, 

’ 

for  the  special  case  of  a lossless,  nondispersive  delay  line  with  lightly  coupled, 
nonrcflccting  taps  (5).  Since  the  signals  in  the  delay  line,  at  the  tap  outputs,  and 
at  the  output  of  the  summing  circuit  are  time-varying,  they  can  be  filtered  in  the 
time  domain  as  well  as  in  the  space  domain. 
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The  tap  weights,  w^  can  be  regarded  as  the  sample  values,  w(r^) , of  a 
weighting  function,  w(t) , where  tj  is  the  delay  to  the  i tap.  These  weights  can 
also  be  regarded  as  the  sample  values,  v(£^),  of  a corresponding  space  weighting 
function,  v(£)  = w(£/c),  where  £.  is  the  distance  to  the  i**1  tap  and  c is  the  propa- 
gation velocity  of  waves  on  the  line.  It  is  assumed  that  all  taps  have  the  same 
impulse  response,  u(t).  Thus,  the  time-varying  function  that  appears  at  the  out- 
put of  the  i**1  tap  when  an  impulse  is  applied  to  the  line  is  u(t)  weighted  by  w.  and 
delayed  by  t..  The  impulse  response  of  the  filter,  h(t),  is  the  summation  of  these 
individual  impulse  responses. 

h(t)  = £ w jU (t  - T i) • (1) 

For  all  of  the  implementations  of  transversal  filters  to  be  discussed  in  this 
paper,  the  taps  have  a bandpass  filter  characteristic,  and  the  reconstruction  of 
w(t)  from  the  sample  values,  wlt  is  accomplished  by  the  filtering  action  of  the 
taps,  or  by  this  action  combined  with  that  of  a simple  compensating  filter  at  the 
input  or  output  of  the  transversal  filter.  If  u(t)  is  redefined  to  include  this  com- 
pensation, when  required,  then  the  right  hand  side  of  (1)  is  identical  in  form  to 
the  expression  representing  reconstruction  of  the  function  w(t)  from  the  samples 
wt  [2],  and  (1)  reduces  toh(t)  = w(t). 

The  output  of  a transversal  filter  for  a general  input  signal  is  then  the  con- 
volution of  that  signal  with  the  filter's  tap-weighting  function: 

y(t)  = [x*w](t)  (2) 

where  y(t)  is  the  output  signal,  x(t)  is  the  input  signal,  and  * denotes  convolution. 
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Thus  transversal  filter  synthesis,  in  terms  of  its  time-domain  specification, 
h(t),  is  extremely  simple.  The  synthesis  procedure  consists  merclj-  of  sampling 
the  specified  h(t)  at  intervals  correr ponding  to  its  bandwidth,  using  ther  o sample 
values  as  weights  on  appropriately  spaced  delay-line  taps,  and  summing  (he  out- 
puts of  these  weighted  taps. 

Scatterer  Distribution  Mapping 

As  an  example  of  the  applications  of  transversal  filters  to  linear  signal  pro- 
cessing, the  mapping  of  scatterers  distributed  in  range  and  doppler  will  be  de- 
scribed1. Such  a mapping  is  of  interest  in  sonar  [6],  radar  [7],  and  radar  astron- 
omy fS] . A finite-energy  signal,  s(t),  is  transmitted  and  the  scattered  energy  is 
received  in  the  presence  of  white  Gaussian  noise.  For  die  detection  of  s(l)  in 
such  noise,  both  a matched  filter  and  a correlator  are  optimum  [9],  A matched 
filter  (Fig.  2)  is  a linear  filter  whose  impulse  response,  h(t),  is  the  time-reversed 
replica  of  the  signal  to  be  detected:  h (t)  = s(-t).  A correlator  (Fig.  3)  multiplies 
the  input,  consisting  of  signal  and  noise,  by  a locally  formed  replica  ol  the  signal 
and  integrates  the  product. 

The  mapping  of  a distribution  of  scatterers  is  a more  complex  problem  than 
that  of  simple  signal  detection;  in  this  case,  the  matched  filter  and  correlator 
are  not  always  optimum. 


1The  echo  location  in  range  and  Doppler  of  a point  or  distributed  targe*  is  a 
special  ease  of  scatterer  distribution  mapping. 
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Nevertheless,  even  when  they  are  not  optimum,  they  often  provide  acceptable 
performance  and  when  they  do  not,  slightly  mismatched  versions  may  provide 
nearly  optimum  performance  [10].  As  a consequence,  these  processors  are  used 
extensively. 

The  mapping  of  a distribution  of  scatterers  by  either  of  these  processors  can 
be  expressed  in  terms  of  the  resolution  function,  \(t,v).  This  function  rep  -esents 

the  response  of  the  processor  when  matched  to  s(t)  for  an  input,  s (t),  which 

T,  V 

is  a version  of  s(t)  delayed  by  r,  with  Doppler  v,  where  v = 2(R/c)f0,  R is  the  range- 
rate,  c is  the  propagation  velocity,  and  f0  is  the  mean  frequency.S * * 8  Since  \(t,  v) 
is  the  response  of  a processor  matched  to  s(t)  = s0,0(t),  it  is  the  convolution  of 

ST,i/With  S°'°; 

X(T,u)  = s *s"0  (3) 

1,1/ 


SThe  approximation  made  here,  that  v = 2(lVc)f0),  is  valid  when  the  total  range- 
rate  interval  of  interest  is  small  relative  to  c/2TvV.  In  this  case,  Doppler  can 
be  approximated  as  a constant  frequency  shift  of  all  the  signal's  spectral  com- 
ponents. This  is  the  case  which  is  usually  treated  in  the  literature,  where  the 
uncertainty  function,  |xl,  and  the  ambiguity  function,  |x|8,  arc  studied  rather 
than  x-  When  the  range-rate  interval  is  not  small  relative  to  c/2TW,  the 
Doppler  shift  of  each  spectral  component  is  proportional  to  the  component's 
frequency  and  the  mapping  becomes  more  complicated  than  the  simple  version 

given  here. 

G 

....  . 


where  s“  (1)  = s (-1). 

T > V T,V 

To  map  the  point  t0» v0  in  the  range-Doppler  space,  the  processor  is  matched 
not  to  s(t)  but  to  st0  ,y0  (t).  However,  the  processor  when  so  matched  responds 
not  only  to  the  signal  scattered  from  the  point  t0,  v0,  but  also  to  scattering  from 
other  points  in  the  t-v  space.  The  output  of  the  processor  for  scattering  from  the 
point  r,  v,  when  matched  to  die  point  rc,  v0,  is  \0(r,v)  = s * s"  . v0(r,p)  is  a 
version  of  \(t , V)  which  has  been  translated  from  the  origin  to  the  point  r0,  vc  and 
has  suffered  a slight  distortion  in  shape  in  the  process  of  translation.  If  x0  is 
concentrated  in  the  region  of  T0,  v0,  or  equivalently  if  x is  concentrated  at  the 
origin,  then  the  processor  output  will  be  a good  estimate  of  the  scattering  cross 
section,  avei’aged  over  the  region  of  concentration.  Thus  the  size  of  the  region 
around  the  origin  in  which  x is  concentrated  determines  the  resolution  with  which 
die  t-v  space  is  mapped  [11]  - [13].  The  nominal  resolution  for  a signal  with 
bandwidth  W and  duration  T is  l/W  in  t and  l/T  in  v. 

The  map  of  the  scattering  distribution  that  results  from  such  processing  is  a 
function  of  r and  v and  can  be  regarded  as  having  bandwidth  in  both  these  variables; 
the  nominal  bandwidth  in  t is  W and  in  v is  T.  A straightforward  extension  of  onc- 
dimcnsional  sampling  theory  to  two  dimensions  leads  to  the  result  that  TW  samples 
per  unit  square  in  the  t- v space  (i.  e. , one  sample  for  each  resolution  cell  of  size 
l/W  in  r and  l/T  in  v)  span  the  entire  space.  The  continuous  map  of  the  scatter- 
ing distribution  can  be  recovered  from  these  samples  by  appropriate  smoothing. 

Although  the  nominal  resolution  in  t and  v is  l/W  and  l/T  for  any  signal 
with  bandwidth  W and  duration  T,  the  theory  of  signal  design  leads  to  many  choices 
of  signal  form,  with  corresponding  variations  in  the  form  of  \(t,v).  A selection 
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from  among  these  choices  must  be  made  on  the  basis  of  prior  knowledge  of  the 
scattering  distribution  as  well  as  the  desired  mapping  resolution.  Signal  design 
is  an  active  area  of  current  research  [11].  The  mapping  may  require  any  signal 


with  time-bandwidth  product  TW,  if  the  resolution  in  the  T,  v space  is  to  be  6^  = 
l/W  and  6 = l/T.  If  the  filter  synthesis  procedure  is  to  be  sufficiently  versatile 

to  provide  a processor  with  the  required  resolution  for  the  mapping  of  any  scattering 
distribution  encountered,  then  this  procedure  must  be  capable  of  synthesizing  a 
filter  to  match  any  signal  with  this  time-bandwidth  product.  When  the  TW  pi’oduct 
becomes  large,  the  frequency  domain  synthesis  of  conventional  longitudinal  filters 
matched  to  arbitrary  signals  specified  in  the  time  domain  becomes  increasingly 
difficult.  On  the  other  hand,  the  synthesis  of  transversal  filters  remains  simple, 
since  the  procedure  consists  merely  of  selecting  the  tap  weights  to  cori’cspond  to 
the  sample  values  of  the  impulse  response. 

Suppose  it  is  desired  to  map  a distribution  of  scatterers  in  the  range  interval, 

• • • 

4' It  - RnaX  - Raln,  and  the  range-rate  interval,  Alt  = Raax  - Iisln,  with  a prescribed 
resolution  in  range,  6R,  and  range-rate,  65.  Since  t»R  = cfi^/2  and  6;  = Xo  hj 2, 
the  required  signal  duration  and  bandwidth  are 

T = Xo/26j  and  W = c/26R  (4) 

where  \0  is  the  wavelength  corresponding  to  the  mean  frequency,  f0.  From  the 
second  equation  of  (1),  a relation  can  be  established  between  the  system’s  frac- 
tional bandwidth  and  the  number  of  range  resolution  cells  per  unit  wavelength: 

W/f o - X</26r  (5) 
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where  W/f0is  the  fractional  bandwidth.  From  (4),  the  required  si  , n il 
time-bandwidth  product  is  given  by 

TW  = cXoAl&A  (6) 

The  time  required  to  form  the  receiver  output  for  a given  rang:  , U,  and  range- 
rate,  R,  within  a resolution  cell  of  size  6R  * t)R  is  the  resolution  cell  processing 
time,  tp.  The  duration  of  the  reflection  from  scatterers  within  the  range  interval 
AR  is  (2  Ar/c)  + T,  and  the  effective  resolution-cell  processing  time  for  a matched 
filter,  tv[0f],  is  this  duration  divided  by  the  number  of  range  resolution  cells, 

Ar/Gr,  or  tpjnf  j = [(2A«/c)  + T]  6r/£.'r.  However,  for  the  correlator,  the  resolu- 
tion-cell processing  time  is  simply  the  signal  duration,  tp[eor]  •-  T.  The  time 
between  transmissions  is  the  pulse  interval,  tlf  and  is  available  for  signal  pro- 
cessing,  so  that  the  maximum  number  of  computations  possible  between  trans- 
missions is  Np  = ti/tp,  and  depends  on  whether  a matched  filter  or  a correlator 
is  being  used: 


Np[nf]  = tl/[26R/c)  + (TV  AOl 

(7a) 

Np[eor]  = tj/T. 

(7b) 

The  number  of  computations,  Nr,  required  to  span  an  aroa  £R  * in  the  range 
range-rate  space  is  the  product  of  the  number  of  range  reso;  lion  cell:.,  AR/6K, 
and  the  number  of  range- rate  resolution  cells,  Nr  = (£RA\  ) (/  */6R). 

The  required  number  of  computations  can  be  executed  during  the  interval  tj  only 
if  Np/Nr  > 1. 
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Since  Nr  is  the  same  for  both  the  matched  filter  and  the  correlator,  the  ratio 
Np[af]/Nv[cer]  is  a measure  of  the  processing  speed  of  the  matched  filler  relative 
to  the  correlator.  From  (7a)  and  (7b), 

Np[of]/Np[cor]  ~ & TW  (8) 

where  a = l/[l  + (cT/2Ar)]  and  0 < a < 1.  Only  when  a < l/TW,  does  the  cor- 
relator have  greater  processing  speed  than  Ihe  matched  filter. 

These  ideas  will  now  be  illustrated  with  numerical  examples  from  sonar  and 
radar  astronomy.  Although  the  examples  arc  not  intended  to  correspond  to  actual 
systems,  they  are  representative  of  the  parameters  that  might  be  required  for 
ocean-floor  or  planetary  mapping.3  The  system  parameters,  along  with  the  signal 
parameters  required  to  achieve  these  system  parameters,  arc  presented  in 
Table  I,  and  the  corresponding  processor  performance  parameters  are  presented 
in  Tabic  II.  Note  that  in  the  sonar  example,  AR  is  not  small  l'elative  to  c/2TYV 
and,  therefore,  the  frequency  shift  approximation  for  Doppler  is  not  valid  over 
the  entire  range-rate  interval.  However,  the  nominal  range-rate  resolution  is 
still  l/T.  A method  for  treating  this  case  by  dividing  AR  into  subintervals 
which  are  small  relative  to  c/2TW  is  discussed  later. 


3 Since  both  the  general  expressions  and  the  numerical  values  arc  shown,  other 
systems  can  be  compared  either  by  direct  calculation  or  by  scaling. 
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Table  I 


Symbol 

Sonar 

Radar  Astronomy 

SYSTEM  PARAMETERS 

Wavelength 

Xo 

10  cm 

10  cm 

Propagation  Velocity 

c 

1.5  x 103  m/s 

3 x 108  m/s 

Mean  Frequency 

fo 

15  kHz 

3 GHz 

Range  Resolution 

6r 

1 m 

2.5  * 10*  m 

Range-Rate  Resolution 

6r 

0.05  m/s 

2 x 10"3  m/s 

Range  Interval 

103  m 

107  m 

Range-Rate  Interval 

5 m/s 

2 m/s 

Pulse  Interval 

tl 

3 s 

3 min 

SIGNAL  PARAMETERS 

Signal  Duration 

T 

1 s 

25  s 

Bandwidth 

W 

750  Hz 

6 kHz 

Time- Bandwidth  Product 

TW 

750 

1.5  x 105 

Fractional  Bandwidth 

w/f0 

0.05 

2 x 10-6 

From  Table  I,  the  bandwidth  and  TW  product  for  the  sonar  example  are 
750  Hz  and  750,  and  for  the  radar  astronomy  example  are  6 kHz  and  1.5  x 105, 
respectively.  Ultrasonic  transversal  filters  have  been  constructed  with  bandwidths 
to  1 MHz  and  TW  products  to  103,  and  filters  are  under  development  with  bandwidths 

4 

to  1 GHz  and  TW  products  to  10  [15],  [16],  Tiius,  existing  filters  have  bandwidths 
in  excess  of  the  requirements  of  sonar  and  radars  with  system  parameters  com- 

4 

parable  to  those  assumed  in  the  examples.  'Hie  TW  product  of  these  existing 
filters  is  adequate  for  such  sonars.  Since  the  TW  product  is  a measure  of  the 
maximum  possible  mapping  resolution,  sonar  resolution  could  be  inci'eased  by 
more  than  an  order  of  magnitude  beyond  the  value  assumed  in  this  example  with- 
out exceeding  the  TW  capacity  of  transversal  filters  currently  under  development. 
Conversely,  even  these  latter  filters  have  TW  capacities  which  are  an  order  of 
magnitude  smaller  than  that  required  for  the  radar  astronomy  example.  It  is 
hoped  that  advances  in  the  technology  of  microwave  acoustic  delay  lines  will  pro- 
vide the  increase  in  TW  capacity  needed  for  such  radar  astronomy  applications. 

From  Table  II,  the  speed  with  which  the  matched  filter  can  map  the  distribution 
of  scatterers  is  approximately  400  times  that  of  the  correlator,  for  both  the  sonar 
and  radar  astronomy  examples.  However,  the  values  of  Np/Nr  presented  in 
Table  U indicate  that  even  the  matched  filter  is  approximately  100  times  too  slow 
to  accomplish  the  mapping  in  real-time,  for  either  the  sonar  or  the  radar. 

4It  will  be  shown  later  that  excess  bandwidth  can  be  used  to  increase  processing 


The  simultaneous  use  of  a sufficient  number  of  individual  filters,  each  matched 
to  a different  Doppler  modification  of  the  signal,  will  permit  real-time  processing. 

An  alternative  method  that  uses  a single  filter  and  employs  time  compression  and 
serial  processing  will  now  be  described.  If  desired,  a combination  of  these  techniques 
may  be  used. 

If  the  signal  to  be  processed  is  compressed  in  time  by  the  factor  k,  then  the 
resolution  cell  processing  time,  tp,  is  reduced  by  this  factor  and  the  number  of 
computations,  Np,  during  the  interval  tl  is  increased  by  this  factor.  Thus,  time 
compression  by  the  factor  k permits  a k-fold  increase  in  processing  speed.  For 
the  examples  given  earlier,  real-time  processing  can  be  accomplished  with  a 
single  matched  filter  if  k equals  77  for  the  sonar  or  110  for  the  radar. 

The  bandwidth  of  a time-compressed  signal  is  increased  by  the  factor  k and 
the  implementation  of  a matched  filter  for  a time-compressed  signal  must  have  a 
bandwidth  capability  k times  larger  than  the  bandwidth  of  the  corresponding  real- 
time signal.  If,  by  a particular  implementation,  a transversal  filter  can  be  con- 
structed which  is  capable  of  producing  impulse  ’'esjxmses  with  a maximum  band- 
widthof  WBBX  and  a maximum  duration  of  TBax,  the  filter  can  be  used  for  time- 
compressed  processing  of  any  signal  with  bandwidth  W and  duration  T such  that 
W<VVBM  and  TW  s TBaxWBM;  the  maximum  compression  factor  is  kBBB  = WBBX/W. 

Thus,  bandwidth  capability  in  excess  of  the  bandwidth  of  the  real-time  signal 
makes  possible  the  use  of  time  compression.  For  the  sonar  and  radar  astronomy 
example,  the  excess  bandwidth  of  existing  fillers  makes  possible  values  of  k 
large  enough  to  permit  real-time  processing  by  a single  filter. 
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Since  kEaX  is  proportional  to  WB<X,  then  WBBJ(  is  a figure  of  merit  for  tlie  filter's 
processing  speed,  and  since  the  maximum  possible  number  of  resolution  cells  per 
unit  cell  is  proportional  to  the  signal's  time -bandwidth  product,  TW,  then  TB(U(WBt>x 
is  a figure  of  merit  for  the  filter's  resolution  capabilities.  Using  current  ultra- 
sonic technology,  fillers  and  time  compressor’s  have  been  built  which,  for  most 

sonar  systems,  have  a kBBX  greater  than  103  [17],  and  microwave  acoustic  devices 

6 

arc  under  development  that  will  make  possible  values  of  kBax  greater  than  10  [18]. 
However,  in  many  radar  systems  the  signal  bandwidth  limits  kBiUt  to  small  values 
for  these  same  filters.  It  is  hoped  that  additional  developments  in  microwave 
acoustics  will  soon  make  larger  time  compressions  possible. 

A Generalized  Linear  Receiver 

The  use  of  time  compression  with  matched  filtering  is  illustrated  by  the 
linear  receiver  shown  in  Fig.  4.  The  matched  filter  shown  in  Fig.  4 has  multiple 
output  ports,  whose  significance  will  be  discussed  later.  A range -delayed  and 
Doppler-modified  signal,  sT  j,(t),  is  to  be  px’occsscd  by  the  receiver  whose  band- 
width, Wr,  must  be  wide  enough  to  accomodate  the  possible  Doppler  shifts  of  the 
signal's  spectral  comjxonents.  By  means  of  singe-sideband  (SSB)  modulation,  the 
receiver  shifts  the  spectrum  of  the  input  signal  downward  by  the  amount 

! 

f4  = f0  - VVr/2,  from  the  receiver's  passband  to  its  baseband,  as  shown  in  Fig.  5a. 

'The  factor  1/2  appears  here  because  the  signal  bandwidth  is  two-sided,  so  that  it 
contains  l>oth  positive  and  negative  frequencies,  w'hercas  for  bandpass  receivers  the 
bandwidth,  \Vr,  is  traditionally  single-sided. 
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The  receiver’s  input  is  so  shifted  b.  sc  Ihc  required  sampling  rate  is  often  lower 
for  baseband  signals  than  for  bandpass  espials  [2].  The -lowest  possible  sampling 
rate  is  desirable,  since  the  time  compressor  chosen  for  illustration  is  a sampled- 
data  device  for  which  the  largest  compression  factor,  kBBX,  is  obtained  when  the 
sample  rate  is  lowest.  If  the  sped)  a i of  the  original  signal  is  S(f),  then  the 
spectrum  of  the  Doppler -modified  v rsion  is  S(/?  f),  where  $ = 1 + 2R/e.  It 
follows  that  the  spectrum  of  the  signal  in  the  baseband  is  S(/5f  - fd)  and  the  spectrum 
of  the  time -compressed  baseband  sip  1 is  S([jgf  - fj/k),  as  shgwn  in  Fig.  5b.  A 

second  SSB  modulator  shifts  the  tim . compressed  signal  upward  by  the  amount 

6 

fu  to  the  frequency  band  of  the  matched  filter.  The  shifted  spectrum  is 

S([(/8f  - fd)/k]  + fu).  'Die  eompressh  i factor  k has  boon  chosen  so  that  after  time 

compression  the  signal  bandwidth  c responds  to  the  filter  bandwidth,  k = k0„. 

If  a single  sample  is  t alien  from  each  cell  of  size  x 6 j,  in  the  range  range-rate 
space,  the  resulting  set  of  samples  is  sufficient  for  an  approximate  description  of 
the  scattering  distribution.  Since  a m {died  filter  processes  all  values  of  range  for 

a single  value  of  Doppler,  this  space  ecu  be  mapped  with  filters  matched  to  the 

• • 

values  of  R at  the  centers  of  strips  of  width  6*  parallel  to  the  R = 0 axis. 

Therefore,  a total  of  A 'H  /6j  individual  matched  filters  are  required  for  l’cal-time 
processing.  When  a filter  operates  in  compressed  time  with  a compression  factor 
k,  compressed  blocks  of  signal  can  be  processed  serially  at  a k-fold  greater  rate. 
Tims,  whenk  serial  processing  can  be  substituted  for  parallel  processing 

and  only  a single  filter  is  required. 

6 This  shift  is  not  required  if  the  mat'  bed  filter  operates  at  baseband. 


Time  compression  and' transversal  filler :ng  are  the  features  that  distinguish 


this  receiver  from  a conventional  superhetrodyne  receiver  which  employs  only 
frequency  translation  and  bandpass  filtering.  Frequency  translation,  usually 
accomplished  by  SSB  modulation,  changes  the  fractional  bandwidth  of  the  input 
while  leaving  its  duration  unclm  whereas  time  compression  changes  the 
duration  while  leaving  the  fractional  bandwidth  unchanged.  Thus  this  generalized 
receiver  provides  complete  flc'ib'lily  in  utilizing;  the  performance  characteristics 
of  the  matched  filter  and  mala  p >ssiblc  the  efficient  application  of  transversal 
filters,  implemented  with  mic:t:o>.n  re  acoustic  devic  es,  to  the  processing  of  signals 
with  spectra  anywhere  in  the  fr.  qnency  range  front  seismic  through  radio  to  micro- 
wave  frequencies. 

Time  Compression 

Time  compression  is  achieved  by  storing  the  signal  in  memory  and  retrieving 
it  at  a rate  greater  than  the  rate  of  storage  [19],  Although  time  compression  can  be 
achieved  with  analog  devices  operating  on  continuous  signals  [20],  [21],  sampled  data 
devices  arc  usually  employed.  The  principle  of  time  compression  can  be  illustrated 
by  a transversal  filter  with  p uniformly -weighted  taps  (Fig.  6)  which  are  separated 
by  time  intervals  of  duration  /'  rt  - At , - Arc , where  Arr  is  the  sample  interval 
in  real-tune,  Atc  is  the  sample  interval  in  compressed  time,  and  p is  an  integer 
greater  than  or  equal  to  T\V.  A!  the  instant  when  the  oldest  sample  in  the  line  passes 
the  Iasi  tap,  it  appears  at  the  output  of  the  filter.  Atc  seconds  later,  the  next-to- 
oldcst  sample  passes  the  next  to  last  tap  and  appears  at  the  output  of  the  filter, 
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Fig.  G.  Transversal  Filter  lime  Compressor,  p = 4.  Relative  positions 
in  space  of  the  output  samples  are  shown  on  the  right  for  specified 
time  delays:  (a)  t = 0;  (b)  t = Atc;  (c)  t = 2Atc;  (d)  t = 3At.  ; 

(e)  t = Arr;  (f)  t = 2Arr  + 3Atc  . 


and  samples  continue  to  pass  their  corresponding  taps  and  appear  at  the  output 
until  (p  - 1)  At,  seconds  later,  at  which  time  the  newest  sample  passes  the  first 


tap  and  appears  at  the  output  of  the  filter.  Thus  p samples  have  appeared  at  the 
output  of  the  filter  in  the  proper  sequence,  but  with  a separation  of  A tc  instead  of 
At,.;  they  have  been  compressed  in  time  with  a compression  factor  k = Att/Atc. 

A pause  of  duration  (k-p+1)  Atc  follows,  during  which  the  next-to-oldest  sample 
propagates  the  remaining  distance  to  the  last  tap.  The  oldest  of  the  p samples 
has  passed  out  of  the  delay  line,  all  samples  have  propagated  a distance  corres- 
ponding to  A rr  along  the  line,  and  a new  sample  has  been  introduced  into  the 
line.  The  process  just  described  now  repeats.  Thus  the  output  of  the  time 
compressor  consists  of  blocks  of  p compressed  samples  v.  ith  a gap  of  duration 
(k-p+l)ATc  between  blocks.  Each  block  contains  all  of  the  samples  from  the 
previous  block  except  the  oldest,  with  the  samples  shifted  within  the  block  by 
Atc  relative  to  the  positions  of  the  samples  within  the  previous  block,  and  with 
a new  sample  added.  These  blocks  can  be  reconstructed  to  form  blocks  of 
continuous  signal  by  the  filtering  action  of  the  taps  themselves,  if  the  tap  im- 
pulse response  is  properly  selected.  The  description  of  the  operation  given 
here  is  valid  only  when 

(l-l/p)ATr  <Art  < ATr  (9) 

When  this  inequality  is  not  satisfied,  the  samples  do  not  pass  their  corresponding 
taps  sequentially  in  time  and  the  situation  becomes  complicated. 

The  compressed  signal  lies  entirely  within  one  or  more  of  these  blocks  only 
when  the  input  signal  lies  entirely  within  the  delay  line.  Therefore  if  p^TW,  each 
of  q adjacent  blocks,  where  q - p-TWU,  will  contain  the  entire  signal;  all  other 
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blocks  will  contain  only  part  of  the  signal  or  none  at  all.  A filter  matched  to  a 
particular  Doppler  modification  of  the  signal  is  connected  to  the  output  of  the  time 
compressor.  At  the  end  of  a time  interval  (p-1)Atc,  equal  to  the  duration  of  one 
of  the  output  blocks,  the  impulse  response  of  the  filer  can  be  altered  to  match  a 
different  Doppler  modification  of  the  signal.  After  another  At,,  seconds  have  elapsed, 
the  filter  can  again  be  modified  to  match  another  Doppler  modification.  This  process  is 
repeated  at  intervals  of  At,,  until  all  required  Dopplers  have  been  examined  by  a cor- 
responding matched  filter.  Since  there  ai'o  only  q blocks  that  contain  the  entire 
signal,  at  most  q modifications  of  the  filter's  response  can  be  made  before  it  must 
be  returned  to  its  initial  state.  Otherwise,  the  filter  response  that  matches  the 
signal  may  not  be  formed  during  the  time  that  the  signal  is  present  in  the  time- 
compressor  output. 

The  time-compression  factor,  k,  acts  as  a fundamental  limit  on  the  number  of 
Dopplers  that  can  be  searched  for  all  range  values.  The  utilization  factor,  qAi 
measures  the  degree  to  which  a time-compressor's  performance  approaches  this 
fundamental  limit.  In  the  compressor  used  for  illustration,  k must  be  greater  than 
or  equal  to  p.  Otherwise,  the  sample  blocks  would  overlap  and  fewer  than  q Dopplers 
could  be  searched.  The  optimum  situation  exists  when  p can  be  made  equal  to  k, 
for  in  this  case  the  gap  between  blocks  reduces  to  Atc  and  q has  the  maximum  value 
possible  for  a fixed  value  of  k.  ’then  q/k  = 1 - (TW-1  )/k  and,  as  k increases,  the 
utilization  factor  approaches  unity.  Thus  for  the  time-compressor  illustrated  here, 
both  large  k and  large  p arc  desirable. 
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For  most  sonar  systems,  values  of  k greater  than  TW  are  necessary'  if  all 
required  Dopplers  arc  to  be  searched  for  all  ranges,  and  time  compressors  with 
sufficiently  large  k have  been  developed  for  such  systems  [19].  However,  for  many 
radar  systems,  the  required  time -compression  factors  are  less  than  TW  and  the 
time  compressor  choosen  here  for  illustration  is  not  appropriate.  For  such  radars, 
and  for  sonars  with  k<TW,  time  compressors  with  tap  intervals  that  do  not  satisfy 
(9)  must  be  used  or  other  equivalent  devices  must  be  employed. 

It  may  not  be  necessary  to  alter  the  filter's  response  q times  to  search  q 
Dopplers.  Within  a Doppler  interval  that  is  small  relative  to  c/2TW,  the  modi- 
fication of  the  signal  as  a result  of  Doppler  can  be  approximated  by  a frequency' 
translation  of  the  signal.  If  the  Doppler  interval,  Zs*  , satisfies  this  condition,  the 
search  in  Doppler  can  be  accomplished  by  shifting  t he  frequency  of  the  variable - 
frequency  oscillator  (VFO)  of  the  second  SSB  modulator  by'  amounts  corresponding 
to  6^  , at  time  intervals  £rr,  instead  of  by  altering  the  filter's  response.  It  is  still 
true,  of  course,  that  only  q Dopplers  can  be  examined.  If  Aj  is  loo  large  for  the 
frequency  shift  approximation  to  be  valid  over  the  entire  interval,  it  can  be  divided 
into  subintervals  in  which  the  frequency -shift  approximation  is  valid.  Then  the  filter 
must  be  altered  to  match  the  Doppler  modification  of  the  signal  for  the  Doppler  value 
at  the  center  of  each  of  these  subintervals;  within  the  subintervals,  however,  the 
Doppler  search  can  be  achieved  by  frequency  shifts  of  the  VFO. 

Although  this  description  of  time-compressed  processing  is  based  on  the 
assumption  that  the  q different  modifications  of  the  signal  arc  Doppler  modifications, 
the  generalized  receiver  described  here  is  not  so  restricted.  The  q possible 
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modifications  of  the  filter's  impulse  response  could  also  represent  signal 


modifications  such  as  those  resulting  from  acceleration.  Even  more  generally, 
they  could  represent  any  set  of  q orthogonal  signals,  such  as  those  corresponding 
to  q symbols  of  an  alphabet  for  use  in  communication  or  q words  in  a coded  message. 

Although  any  type  of  memory  can  be  used  for  a time  compressor,  computer- 
type  memories  are  usually  employed;  the  signal  is  sampled  at  baseband,  digitized 
to  form  a binary  number,  and  stored  [19].  The  ultrasonic  delay  line  time  com- 
pressor (DELT1C),  introduced  by  V.  Anderson  in  1956  [22]  for  the  compression  of 
sampled  signals  digitized  at  two  levels,  employed  a computer-type  serial  quartz 
delay-line  memory,  which  was  a forerunner  of  modern  recirculating  memory  stores 
[23].  Configurations  of  the  CELTIC  exist  that  produce  blocks  of  samples  identical 
in  form  and  compression  factor  to  those  produced  by  the  transversal  filter  time 
compressor  (TFTC),  both  for  the  simple  case  illustrated  here  and  for  the  more 
complicated  cases  that  occur  when  (9)  is  not  satisfied.  Thus  there  is  no  difference 
between  the  performance  of  a CELTIC  and  a corresponding  TFTC.  However,  the 
total  delay  tunc  required  in  the  CELTIC  is  approximately  £ Tr,  whereas  in  the 
TFTC  it  is  p£rr.  Note  that  p£rr  is  greater  than  or  equal  to  T,  the  signal's  duration, 
and  the  requirement  that  the  TFTC  have  this  much  delay  makes  it  impractical  fox- 
many  systems.  It  was  introduced  here  for  illustrative  purposes  because  of  its 
simplicity.  Convci-scly,  the  CELTIC  is  practical  for  application  to  many  systems 
because  of  its  relatively  short  delay-time  requirements.  Its  delay-time  advantage 
is  obtained,  however,  at  the  expense  of  requiring  binary  digital  signals  at  the  input, 
since  the  recirculating  delay-line  amplitude  stability  is  generally  insufficient  to 
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preserve  the  sample  amplitudes  of  analog  or  multilevel  digital  signals  during  the 
number  of  recirculations  required. 

To  use  the  DELT1C  for  the  compression  of  analog  or  multilevel  digital  signals, 
the  samples  must  first  be  converted  to  binary  coded  numbers,  and  then  compressed 
either  by  an  appropriate  number  of  DELTICS  operating  in  parallel,  or  by  a single 
DEL  TIC  opcx’ating  serially  at  a correspondingly  inci'cased  spccd[2-1]. 

Multiplc-Output-Porl  Transversal  Filters 

For  many  implementations  of  transversal  filters,  it  is  possible  to  realize  a 
number  of  simultaneous  outputs  with  a single  delay  line  [25].  These  multiple  out- 
puts are  foi'mcd  from  an  array  of  tap  weights,  as  shown  in  Fig.  7.  Multiple-output  - 
poi't  filters  may  provide  savings  in  system  size  and  complexity,  and  some  technique:, 
for  their  implementation  ai'e  discussed  in  the  second  section  of  this  paper. 

Such  a filter  can  l'cplace  a bank  of  filters  operating  in  parallel  to  provide  real- 
time pi'oeessing.  In  addition,  multiple-output-port  matched  filters  ai-e  useful  in 
increasing  the  effective  speed  of  sei'ial  processing  in  compressed  time.  If  such  a 
filter  operates  with  a time -compression  factor,  k,  and  has  m output-ports,  and  each 
port  rcpi'cscnts  a different  one  of  the  impulse  responses  needed  for  the  search  in 
Doppler,  then  the  filler  has  an  effective  processing  speed  of  k x m.  For  a system 
with  fixed  Doppler  search  requirements,  the  addition  of  multiple  output -ports 
reduces  the  time  compression  required  for  Doppler  seai'ch,  and  thus  reduces  the 
bandwidth  capability  required  of  the  filter. 

All  of  the  earlier  statements  concerning  the  application  of  sampling  in  space 
and  time  to  transversal  fillers  apply  also  to  multiple -output-port  transversal 
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filters.  The  tap  weights,  \vy,  can  be  regarded  as  the  sample  values,  w^Tj, 
of  a weighting  fiuiction,  \vi  (t) , where  r:  is  the  delay  to  the  i tap  and  the  sub- 
script j indicates  the  port  at  which  the  output  is  observed.  These  weights  can 
also  be  regarded  as  the  sample  values,  Vj(4j),  of  a corresponding  space  weight- 
ing function,  Vj(£)  = Wj(£/c),  where  £1  is  the  distance  to  the  i*  tap  and  c is  the 

til 

propagation  velocity  of  waves  on  the  line.  The  j output  has  impulse  response 

h j(t)  = Lw„u(t  - Tj)  (10) 

1 

where  \vl}  is  the  weight  of  the  i^  tap  for  the  j^1  output  and  u(t)  is  the  tap  impulse 

response.  The  filtering  action  of  the  taps,  with  compensation  as  required,  reconstructs 

Wj(t)  from  the  sample  values  wtJ,  and  (10)  reduces  to  hj(t)  = Wj(t).  Thus,  for  a 

til 

general  input  signal,  the  j output  of  a multiple -output-port  transversal  filter  is 

til 

the  convolution  of  that  signal  with  the  filter's  j tap  weighting  function: 


y j(t)  = [x  *wj  (t) 


(ii) 


.th 


where  y^t)  is  the  j output  signal  and  x(t)  is  the  input  signal.  The  synthesis  of 
multiple -output  -port  transversal  filters  is  comparable  in  simplicity  to  that  for  the 
single -output-port  filter. 

Conclusion 


Because  of  the  extreme  simplicity  of  its  synthesis,  the  transversal  filter  is  ideal 
for  applications,  such  as  linear  signal  processing,  where  a complicated  impulse  re- 
sponse with  large  a time-bandwidth  product  may  be  required.  Existing  ultrasonic 
delay  line  transversal  filters  have  many  applications  in  sonar,  radar  and  radar  astronomy 
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systems,  and  microwave  ultrasonic  delay  line  transversal  fillers  under  development 
can  be  expected  to  have  many  additional  applications.  Since  generally  the  physical 
size  of  these  devices  decreases  and  the  processing  speed  increases  as  the  maximum 
bandwidth  increases,  it  is  to  be  hoped  that  the  operation  of  ultrasonic  delay  line 
transversal  filters  will  be  advanced  still  further  into  the  microwave  region. 

In  many  systems,  real-time  signal  processing  is  impractical  because  the 
number  of  simultaneous  operations  to  be  performed  exceeds  the  system's  capability. 
If  the  bandwidth  capability  of  the  signal-processing  filter  exceeds  the  bandwidth 
of  the  signal  to  be  processed,  the  signal  can  be  time-compressed  by  an  amount, 
k,  proportional  to  the  bandwidth  ratio.  When  the  input  signal  is  thus  compressed 
in  time,  the  processing  can  proceed  at  a k-fold  greater  rate.  Thus  linear  filters 
operating  in  compressed  time  make  possible  significant  increases  in  the  effective 
speed  with  which  linear  signal  processing  can  be  performed,  and  transversal 
filters  arc  ideal  for  this  application  because  of  the  simplicity  of  the  synthesis 
procedure  required  to  obtain  an  arbitrary  impulse  response. 

Some  implementations,  to  be  described  in  the  second  section  of  this  paper, 
make  possible  the  construction  of  transversal  filters  with  multiple  output -ports, 
at  each  of  which  appears  simultaneously  the  response  to  an  impulse  at  the  input- 
port,  each  response  corresponding  to  one  of  a set  of  mutually  orthogonal  functions. 

If  the  filter  has  m such  ports,  and  processes  at  a k-fold  increased  rate,  as  in  the 
generalized  receiver  described  here,  it  is  equivalent  to  k x m single-output -port 
filters  operating  in  parallel  in  real-time. 
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It  should  not  be  inferred  that  the  uses  of  transversal  filters  discussed  so  far 


are  the  only  ones  for  which  they  are  suited.  The  case  with  which  time-domain 
synthesis  of  transversal  filters  with  large  T\V  product  can  be  accomplished,  and 
the  availability  of  compact  implementations  employing  microwave  acoustic  delay 
lines,  make  transversal  fillers,  in  conjunction  with  time  compression  and  frequency 
translation,  useful  in  a wide  range  of  applications  encountered  in  diverse  disciplines. 
Besides  the  many  applications  to  scatterer-distribution  mapping  in  the  fields  of  radar, 
radar  astronomy,  and  sonar,  there  are  also  mapping  applications  in  seismology, 
nondestructive  testing,  and  medical  instrumentation.  In  addition,  there  are  signif- 
icant applications  in  interferometry,  spectral  analysis,  communications,  adaptive 
filtering,  and  delay-line  technology,  as  well  as  in  laboratory  instrumentation. 

A number  of  implementations  of  transversal  filters  and  time  compressors 
employing  ultrasonic  delay  lines  arc  currently  available,  some  of  which  are 
discussed  in  the  second  section  of  this  paper.  Although  these  devices  have  made 
possible  the  solution  of  a number  of  important  signal -processing  problems,  there 
are  still  many  important  problems  for  which  ultrasonic  transversal  fillers  and  time 
compressors  with  higher  speed  (W)  and  larger  capacity  (T\V)  are  required.  Because 
of  the  frequencies  at  which  these  needed  devices  must  operate,  they  fall  within  the 
domain  of  microwave  science  and  engineering.  The  special  talents  of  the  worker  in 
microwaves  might  fruitfully  be  applied  in  advancing  this  technology’. 
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ULTRASONIC  TRANSVERSAL  FILTERS 


In  the  first  section  of  this  paper,  the  Kallmann  transversal  filter  was  de- 
scribed and  some  of  its  properties  and  applications  were  discussed.  In  the  sec- 
ond section,  a generalized  transversal  filter  is  introduced  and  analyzed,  some 
features  of  delay  line  transduction  and  propagation  are  discussed  in  the  context 
of  transversal  filter  implementation,  and  a number  of  practical  implementations 
employing  magnetostrictive,  piezoelectric,  and  optoacoustic  taps  on  a variety  of 
ultrasonic  delay  lines  arc  described. 

The  G ( moralize d Transversal  Filter 

In  the  Kallmann  filler,  it  is  assumed  that  the  desired  signal  is  propagating 
along  the  delay  line  where  taps,  which  do  not  perturb  the  propagating  wave,  spa- 
tially sample,  weight,  and  sum  the  signal.  No  consideration  was  given  to  the 
methods  of  transduction  by  which  the  wave  was  launched.  There  is  often  an 
advantage  in  using  a transducer  composed  of  a weighted  array  of  elements.  If 
its  elements  are  reciprocal,  the  array  can  be  used  both  for  launching  and  tapping, 
which  is  a useful  feature  in  some  applications.  The  utility  of  filters  with  multi- 
ple output  ports  was  discussed  in  the  first  section  of  the  paper,  where  a method 
employing  a two  dimensional  weighting  function  applied  to  a tap  array  was  de- 
scribed. Another  method  uses. a sequence  of  weighted  arrays  along  the  line.  A 
natural  extension  ;s  a filler  with  both  multiple  output  ports  and  multiple  input 
ports.  'Ihe  most  general  filter  would  be  composed  of  a sequence  of  nonoverlapping 
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arrays  along  a delay  line  with  absorbers  at  the  ends  of  the  line  to  prevent  reflec- 
tions, with  a two  dimensional  weighting  function  applied  at  each  array  location 
to  form  multiple  ports,  and  with  each  array  reciprocal  and  usable  as  both  an 
input  and  an  output  array.  Since  several  of  the  implementations  of  transversal 
filters  described  here  have  one  or  more  of  these  features,  the  generalized  trans- 
versal filter  incorporating  all  of  them  is  discussed.  First,  however,  an  applica- 
tion of  such  a filter  with  three  sequential  ports  is  described. 

A filter  that  is  matched  to  a time-reversed  version  of  the  transmitted  signal 
is  frequently  used  for  detection  or  scatterer  distribution  mapping.  The  genera- 
tion of  the  signal  to  be  transmitted  is  straightforward  when  the  signal  is  simple, 
but  when  it  is  complicated  and  has  a large  TW  product,  sophisticated  techniques 
of  generation  may  be  required.  One  method  is  to  sample  the  signal  and  store 
the  samples  in  memory,  to  be  retrieved  and  reconstructed  into  the  signal  at  each 
transmission.  Another  method  is  to  generate  the  signal  recursively  by  the  appli- 
cation of  an  appropriate  algorithm  to  an  elemental  signal.  For  signals  with 
large  TW  product,  the  former  method  requires  a large  memory,  as  does  the 
latter  method  because  of  its  complicated  algorithm. 

An  attractive  alternative  is  to  use  the  same  generalized  transversal  filter 
for  both  signal  generation  and  signal  processing.  Consider  a transversal  fil- 
ter with  two  transducers  separated  by  a tap  array.  If  h(t)  is  the  response  at 
the  output  of  the  tap  array  when  an  impulse  is  applied  to  one  of  the  transducers, 
then  h(-t)  is  the  response  when  the  impulse  is  applied  to  the  other  transducer. 
Thus,  if  h(t)  is  the  required  signal  for  transmission,  it  can  be  obtained  at  each 
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transmission  merely  by  applying  an  impulse  to  the  appropriate  transducer;  the 
other  transducer  is  the  input  of  a filter  matched  to  h(t).  Such  a three-port  trans- 
versal filter  is  illustrated  in  Fig.  8.  If  the  transducers  and  taps  are  reciprocal, 
an  alternate  method  is  to  use  the  tap  array  as  the  filter  input,  with  the  trans- 
ducers acting  as  outputs.  If  the  response  at  one  transducer  due  to  an  impulse 
applied  to  the  tap  array  is  h(t),  then  the  response  at  the  other  transducer  is  h(— t), 
and  again  the  filter  can  be  used  for  both  the  generation  and  the  matched  filtering 
of  h(t).  For  the  case  of  reciprocal  transducers  and  taps,  the  impulse  response 
between  any  pair  of  ports  is  independent  of  which  port  of  the  pair  is  used  as  the 
input.  If  the  taps  and  transducers  are  lightly  coupled  and  nonreflecting,  the  im- 
pulse response  between  the  transducers  and  the  tap  array  is  related  to  the  tap 
weighting  function  in  a manner  analogous  to  the  relation  for  the  Kallmann  filter: 
h13(t)  = h3l(t)  w(t)  and  Iv^t)  = h5a(t)  = w(-t),  where  the  subscripts  1,2,  and  3 
indicate  the  left  end  port,  the  right  end  port,  and  the  tap  array  port,  l'espec lively 
(see  Fig.  8),  and  where  w(t)  is  the  tap  weighting  function  expressed  in  the  time 
domain. 

An  additional  advantage  of  using  the  same  filter  for  both  signal  generation 
and  detection  or  mapping  is  that  the  constraints  on  variations  in  the  total  delay 
time  of  the  line  can  be  relaxed.  When  an  independent  signal  generator  is  used, 
the  total  delay  time  of  the  transversal  filter's  delay  line  must  be  held  constant 
to  a small  part  of  the  delay  time  between  taps;  otherwise,  as  the  total  delay  time 
is  changed  by  temperature  variations  or  other  causes,  the  transversal  filter  will 
no  longer  match  the  transmitted  signal.  To  achieve  the  necessary  stability, 
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isopaustic  (constant  delay)  materials  may  be  required.  However,  if  the  same 
transversal  filter  is  used  both  for  signal  generation  and  detection,  then  as  the 
total  delay  time  drifts,  the  matched  filter  tracks  the  generated  signal.  The  only 
constraint  on  the  total  delay  time  drift,  is  that  it  be  small  during  the  time  interval 
between  transmission  and  reception  of  the  echo  or  scattered  signal  when  it  is 
compared  with  the  delay  time  between  taps. 

Before  proceeding  with  the  analysis  of  the  generalized  transversal  filter, 
some  aspects  of  transduction  and  tapping  must  be  examined.  In  the  Kallmann 
filter,  discrete  distributions  of  point  taps  are  weighted  and  summed  to  form  the 
output.  In  the  generalized  transversal  filter,  it  is  desirable  to  consider  continu- 
ous as  well  as  discrete  distributions  of  elements  both  for  launching  and  for  tap- 
ping. 

From  the  theory  of  sampling,  any  continuously  distributed  array  of  elements 
can  be  represented  as  a discrete  arrav.  with  the  element  spacing  approximately 
equal  to  c/VV,  where  c is  the  propagation  velocity  of  waves  on  the  line  and  W is 
the  maximum  bandwidth  of  signals  to  be  propagated  along  the  line.  Conversely, 
a discrete  array  of  elements  with  such  spacing  and  appropriate  reconstruction 
filtering  as  is  provided  by  properly  compensated  elements  is  equivalent  to  a con- 
tinuously distributed  array  with  the  same  weighting  function.  When  elements 
have  this  spacing,  they  are  said  to  be  "contiguous.  " The  element  weights  can  be 
regarded  as  components  of  a vector.  This  vector  is  the  sample  function,  cor- 
responding to  a continuous  weighting  function. 
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The  properties  of  an  individual  element  of  a discrete  array  can  I ’ scribed 
in  terms  of  its  effective  length.  If  an  impulse  is  applied  to  such  an  1 t,  the 
resulting  wave  propagating  along  the  delay  line  is  the  element's  imp  1 response; 
the  length  cr  tliis  wave  is  the  effective  length  of  the  element.  When  . i-nnul.se 
wave  propagating  along  the  line  passes  such  an  element,  the  result'  output  is 
the  element’s  impulse  response;  the  effective  element  length  is  c fir  the  im- 
pulse response  duration.  If  the  element  is  reciprocal  so  that  it  can  ! r -ed  both 
as  a launch  and  as  a receive  element,  it  has  the  same  effective  lenglh  in  both 
modes.  An  element  whose  effective  length  is  substantially  less  thru  c/W  can  be 
approximated  by  a point  element  with  an  impulse  response  that  is  a 0 -h".  function. 
If  the  effective  length  is  comparable  to  c/W,  the  element  can  be  rep.  . rented  as  a 
point  element  whose  impulse  response  is  the  same  as  that  of  the  act’  I element. 

If  the  effective  length  is  substantially  greater  than  c/W,  the  element  c.  a be  repre- 
sented as  an  extended,  distributed  array. 

The  relation  between  the  launch  and  receive  array  weighting  functions  and 
the  generalized  transversal  filter's  impulse  response  is  obtained  by  ; / analysis 
paralleling  that  given  in  the  first  section  of  this  paper  for  the  Kallmann  filter. 

The  analysis  is  developed  for  distributed  launch  arrays  as  well  as  for  distributed 
receive  arrays.  As  in  the  case  of  the  Kallmann  filter,  the  delay  lb  e is  assumed 
to  be  lossless  and  nondispersive  and  the  array  elements  arc  assumed  to  be 
lightly  coupled  and  nonrcflccting.  The  impulse  response  of  (lie  fill*  i is  derived 
by  the  integration  of  continuously  distributed  weighting  functions.  Ib  .ever,  by 
regarding  the  resulting  expressions  cither  as  Sficltjcs  integrals  or  integrals 
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of  generalized  functions,  it  can  be  shown  that  the  expressions  apply  to  discrete 
as  well  as  to  continuously  distributed  arrays. 

Consider  a particular  continuously  distributed  array,  say  the  k array  on 
the  line.  Associated  with  this  array  is  the  weighting  function,  wJk(t).  The  sub- 
scripts indicate  that  wJk(t)  is  the  weighting  function  for  the  jth  member  of  the  set 
of  multiple  ports  associated  with  the  kth  array.  An  impulse  applied  to  a differ- 
ential element  of  this  array  generates  a wave  on  the  line  with  amplitude 
wJk(£/c)uk(t  - [z  - £]/c)d£,  where  £ is  the  coordinate  of  the  element,  uk(t)  is  the 
impulse  response  or  the  element,  z is  the  coordinate  at  which  the  amplitude  is 
measured,  and  c is  the  propagation  velocity  of  waves  on  the  line.  The  element 
generates  waves  propagating  in  both  directions  along  the  delay  line  and  the  direc- 
tion of  increase  for  £ and  z is  chosen  according  to  the  direction  of  propagation  of 
whichever  of  these  two  waves  is  being  considered.  It  is  assumed  that  all  ele- 
ments of  the  array  have  the  same  impulse  response,  uk(t).  The  amplitude  of  the 
wave  at  the  coordinate  z due  to  the  entire  distribution  of  elements  in  the  launch 
array  i s 


a(t,z)  = f wJk(£/c)  uk(t  - [z  - 4]/c)d4  . 
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Next,  suppose  that  z is  the  coordinate  of  a differential  element  of  the  mth 
array.  Associated  with  this  array  is  the  weighting  function,  wlB(l).  The  sub- 
scripts indicate  that  w,a(t)  is  the  weighting  function  for  the  l11’  member  of  the  set 
of  multiple  ports  associated  with  the  m^1  array.  If  an  impulse  wave  passes  this 
element  at  time  zero,  the  response  of  the  weighted  element  is  wl3(z)u„(t)dz, 
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where  uB(t)  is  the  impulse  response  of  the  element.  As  before,  it  is  assumed 
that  all  elements  of  the  m*  array  have  the  same  impulse  response,  uc(t),  not 
necessarily  equal  to  uk(t).  When  the  wave,  a(t,  z),  passes  this  element,  the 
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response  is  w,B(z)  / a(X,z)  ua(t  - X)d\dz.  The  impulse  response  between  the 
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port  of  the  k^1  array  and  the  1th  port  of  the  m^1  array  is  obtained  by  substi- 
tuting (12)  in  this  expression  and  integrating  the  result  over  z: 

hjklm  (t)  - ///:  wJk(£/c)  wlB(z/c)  uk(X  - [z  - £]/c)  u,(t  - X)d£dXdz.  (13) 

The  right-hand  side  of  (13)  is  the  quadruple  convolution, 

hJklB(t)  = [wj*  * wlB  * uk  * u.  | (t) , (14) 

where  wjk(t)  = w Jk(-t)  is  the  time-reversed  version  of  the  weighting  function  for 
the  launch  array  and  * denotes  convolution.  A factor,  cp,  which  is  a consequence 
of  integrating  the  time  domain  weighting  functions  with  respect  to  the  spatial  var- 
iables £ and  z,  has  been  omitted  from  (14)  since  it  can  be  absorbed  in  the  weight- 
ing functions.  This  expression  for  the  impulse  response  is  valid  for  both  contin- 
uous and  discrete  weighting  functions.  For  the  discrete  case,  the  convolutions 
with  respect  to  the  weights  are  sums  instead  of  integrals. 

The  frequency  transfer  function  is  given  by  the  Fourier  transform  of  (14). 

HJkl0(f)  = W;k(f)W„(f)Uk(f)lTB(f),  (15) 

where  the  upper-case  symbols  denote  the  Fourier  transforms  of  the  correspond- 
ing lower-case  symbols  and  the  superscript  * denotes  complex  conjugation.  Thus, 
the  response  between  the  ports  of  a pair  of  arrays  of  a generalized  transversal 
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filter  is  equivalent  to  the  response  of  four  cascaded  filters  which  correspond  to 
the  launch  array  reversed  in  space  or  conjugated  in  frequency,  the  receive  array, 
a launch  array  element,  and  a receive  array  element. 

Ujj  and  UB  set  the  maximum  possible  bandwidth  of  signals  to  be  filtered.  If 
consideration  is  restricted  to  signals  whose  spectral  components  lie  within  the 
limits  imposed  by  Uk  and  Ua,  the  filtering  effect  of  the  array  elements  can  be  ne- 
glected. For  discrete  arrays,  Uk  and  UB  also  set  the  spacing  required  for  contig- 
uous elements.  If  the  elements  arc  so  spaced  and  if  consideration  is  restricted 
to  signals  with  corresponding  bandwidlh  or  less,  then  the  effect  of  the  elements 
can  again  be  neglected.  Thus  for  both  contiguously  tupped  arrays  and  foi  contin- 
uous arrays,  when  the  preceding  conditions  arc  satisfied,  (11)  and  (15)  reduce  to 


hjkla(t)  = |w£*wla|(t)  (16) 

and 

Hjkintf)  - W*(f)Wla(f).  (17) 

From  cither  (16)  or  (17),  it  can  be  seen  that  the  impulse  response,  hJklB(t),  is 
the  cross-correlation  between  the  receive  array  weighting  function,  w,.(t),  and 
the  launch  array  weighting  function,  wJk(t).  As  in  the  case  of  the  Kallmann 
filter,  if  uk  and  u„  serve  as  reconstruction  filters,  the  discrete  distributions  of 
weights  can  be  replaced  by  the  continuous  functions  for  which  the  weights  arc 
samples. 
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For  a general  signal,  x(t),  applied  to  the  j port  of  the  k array,  the  out- 
put. at  the  l^1  port  of  the  m^1  array  is 

y jki*  (t)  = | x * w~  * wlB  I (t) . (18) 

i 1 * 

When  the  elements  of  both  arrays  are  reciprocal,  yjklD  = y,E,k , since  in  that  case 

^Jkla  — ^ln  Jk  • 

The  net  effect  of  the  additional  complexity  of  the  generalized  filter  is  to  cas- 
cade the  filtering  action  of  the  launch  array  with  that  of  the  receive  array.  Thus 
the  analysis  and  synthesis  of  the  generalized  transversal  filter  are  simple  exten- 
sions of  the  analysis  and  synthesis  for  the  Kallmann  filter,  hi  the  synthesis 
procedure,  the  designer  has  a choice  in  the  partitioning  of  a specified  impulse 
response  or  frequency  transfer  function  between  the  launch  and  receive  arrays. 
Some  implementations  which  demonstrate  the  advantage  of  this  flexibility  are 
described  later. 

Propagation  and  Transduction 

The  foregoing  theory  of  transversal  filters  is  based  on  the  assumptions  of 
nondispersivc  propagation  and  lightly  coupled  and  nonreflecting  transducer  and 
tap  elements.  Since  one  of  the  most  attractive  features  of  this  theory'  is  the 
simplicity  of  transversal  filter  analysis  and  synthesis,  these  assumptions  are 
desirable  from  a practical  point  of  view.  When  they  are  not  satisfied,  the  analy- 
sis and  synthesis  become  difficult  (e.g.  , [20]). 

'There  are  a number  of  other  practical  considerations  in  the  implementation 
of  transversal  filters.  The  bandwidth  and  TW-product  requirements  of  the 
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particular  application  are  significant  factors  in  the  choice  of  an  implementation 


for  a transversal  filter.  The  bandwidth  has  an  important  relation  to  the  physical 
size  of  the  filter.  For  most  applications,  the  length  of  an  array  in  the  filter  is 
approximately  equal  to  cT/k,  where  k is  the  time  compression  factor.  Since 
k = Ww/W,  where  W0M  is  the  maximum  possible  bandwidth  for  the  filter,  the 
array  length  becomes  L = cTW/Waax.  Thus,  for  a fixed  TW  product,  the  length 
of  the  array  is  inversely  proportional  to  the  filter's  bandwidth  capability  and 
large  WBM  implies  small  size.  Consequently,  the  filter's  mean  frequency,  f0, 
does  not  affect  its  length;  it  makes  no  difference  whether  a filter  with  a bandwidth 
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of  100  Mllz  is  centered  at  100  MIIz  or  at  1 GIIz.  Since  nearly  all  materials  have 
an  acoustic  attenuation  that  increases  monotonically  with  .frequency,  and  since 
the  complexity  of  the  peripheral  electronics  increases  with  frequency,  small 
values  of  f0  are  desirable.  Small  values  of  fc,  to  achieve  low  loss  and  simple 
electronics,  with  large  values  of  WBax,  to  achieve  small  physical  size,  imply  a 
large  fractional  bandwidth,  WBax/f0.  However,  it  becomes  increasingly  more 
difficult  to  obtain  launch  transduction  with  both  low  insertion  loss  and  low  reflec- 
tivity as  the  fractional  bandwidth  is  increased.  Thus,  a compromise  must  be 
made.  For  most  of  the  ultrasonic  transversal  filter  (UTF)  implementations  to 
be  described  here,  fractional  bandwidths  between  0.5  and  1 are  feasible. 

The  ideal  delay-line  material  for  transversal  filters  has  a dispersionlcss, 
low-loss  mode  of  propagation  and  has  the  capability  of  wide -band,  low-loss  trans- 
duction. The  extent  to  which  these  various  conflicting  requirements  can  be  met 
determines  the  success  of  a particular  implementation.  Sonic  implementations 
which  have  proved  to  be  successful  arc  described  below. 
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In  the  frequency  region  from  1 MHz  to  10  MHz,  satisfactory  transduction, 


1 

propagation,  and  tapping  can  be  obtained  with  a wire  delay  line  propagating  tor- 
sional waves,  using  magnetos trictive  transducers  and  taps;  in  the  frequency 
region  from  10  MHz  to  500  MHz,  they  can  be  obtained  with  a single-crystal  de- 
lay line  propagating  surface  waves,  usually  employing  piezoelectric  transducers 
and  taps;  and  in  the  frequency  region  from  500  MHz  to  2000  MHz,  they  can  be 
obtained  with  a single-crystal  delay  line  propagating  bulk  waves,  either  longi- 
tudinal or  shear,  using  thin  film  piezoelectric  transducers  and  oploacoustic  taps. 

Each  of  these  implementations  arc  discussed  in  detail  later.  In  addition  to 
the  combinations  of  delay-line  materials,  transduction,  propagation,  and  tapping- 
used  for  these  implementations,  there  are  many  others  that  have  been  or  could 
be  used.  Since  propagation,  transduction,  and  tapping  have  been  studied  exten- 
sively for  delay-line  applications,  the  reader  is  referred  to  the  literature  [27]  - 
[30].  Some  important  results  arc  summarized  in  this  section. 

There  are  two  fundamental  waves  in  unbounded,  isotropic  solids:  a dilational 
wave  with  propagation  velocity  cp  - (X  + 2 4)'/2/p'/s»  and  a rotational  wave  with 
propagation  velocity  cs  = pV^/p1/2,  where  X and  p are  the  Lame/ constants  of  the 
material  and  p is  the  density.  The  former  is  usually  called  a longitudinal  wave 
since  the  particle  motion  is  in  the  direction  of  propagation;  the  stresses  in  the 
direction  of  particle  motion  arc  tensions  and  compressions.  The  latter  is 
usually  called  a shear  wave  since  the  particle  motion  is  transverse  to  the  dircc- 
tion  of  propagation,  thus  causing  shear  stresses  in  the  direction  of  particle 


r 
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of  boundaries  complicates  the  situation,  causing  energy  conversion  between  the 
two  fundamental  modes,  with  dispersion  as  the  usual  consequence,  hi  some  cases 
essentially  nondispersive  modes  of  propagation  exist  in  bounded  media.  Of  these 
modes,  the  most  important  for  UTF  implementations  are  (1)  extensional  waves  in 
wires  and  ribbons  whose  dimensions  transverse  to  the  propagation  direction  are 
less  than  a wavelength,  with  propagation  velocity  ce  = pl/2( 3X  + 2 p)^2/p^2(X  + p)1/2 
= E^/p^,  where  E is  Young's  modulus;  (2)  torsional  waves  in  wires  whose 
diameter  is  less  than  a wavelength,  with  c6  = p'/2/p1//p;  (3)  width-shear  waves  in 
sheets  whose  thickness  is  less  than  a wavelength  but  whose  width  is  many  wave- 
lengths, with  c,  = u^s/p^p;  (4)  surface  (Rayleigh)  waves  on  plates  whose  thickness 
is  several  wavelengths,  with  c„  = p’/s(0.87  + 1.12a)/p'/H(l  + a)  where  a = X/2(X  -i  p) 
is  Poisson's  ratio;  and  (5)  longitudinal  or  shear  bulk  waves  in  rods  whose  trans- 
verse dimensions  are  many  wavelengths,  with  velocities  ep  and  c6,  respectively. 

The  attenuation  in  polycrystalline  materials  contains  terms  proportional  to 
the  first  and  proportional  to  the  fourth  power  of  frequency.  The  first-power  term 
represents  the  anelastic  loss  and  for  many  materials  decreases  as  the  elastic 
modulus  increases.  The  fourth-power  term  represents  the  Rayleigh  scattering 
loss  and  increases  with  the  size  of  the  individual  crystallites  [30].  The  scattering 
loss  can  be  reduced  in  some  wires  and  ribbons  by  the  partial  orientation  of  the 
crystallites  resulting  from  cold  working.  In  magnetostrictive  materials,  there 
are  additional  losses  due  to  hysteresis  and  eddy  currents  [31],  and  at  frequencies 
below  10  MHz  these  arc  dominant  losses.  Hysteresis  losses  depend  on  the  mag- 
netic state  of  the  material  and  eddy  current  losses  depend  on  the  magnetostrictive 
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coefficients  and  resistivity.  Magnetostrictive  delay  lines  operating  in  the  tor- 


sional mode,  at  frequencies  of  about  1 MHz,  can  be  constructed  from  wire  with 
moderate  resistivity  such  as  NiSpan  C,  with  losses  as  low  as  a few  tenths  of  a 
decibel  per  meter,  so  that  lines  with  as  much  as  30  meters  of  length  and  a TW 
product  of  104  are  feasible.  An  extremely  simple  and  effective  method  of  tapping 
such  lines  is  described  later. 

A composite  structure  can  be  used  to  reduce  the  magnetic  losses.  A thin 
magnetostrictive  film  is  formed,  deposited,  or  plated  on  a nonmagnetic  base 
wire.  The  base  wire  controls  the  propagation  characteristics  and  the  magneto- 
strictive film  controls  the  tapping  characteristics.  This  type  of  structure  has 
been  used  with  extensional  waves  [32],  [33],  and  appears  promising  for  applica- 
tions in  which  the  impulse  response  of  the  filter  must  be  alterable.  Another 
interesting  composite  structure  consists  of  a thin  magnetostrictive  film  vacuum 
deposited  on  an  amorphous  glass  or  fused  silica  plate  [34].  The  attenuation  in 
these  amorphous  materials  contains  terms  proportional  only  to  the  first  and  sec- 
ond powers  of  frequency,  so  that  the  propagation  losses  arc  acceptable  at  much 
higher  frequencies  than  for  polycrystalline  materials.  Wide-band  delay  lines 
operating  in  the  width-shear  mode  at  frequencies  to  100  MHz  can  be  constructed 
in  this  manner. 

Amorphous  dielectric  plates  can  be  used  also  for  the  propagation  of  surface 
waves  [35].  'this  mode  of  propagation  has  received  much  attention  in  recent  years, 
not  only  for  amorphous  plates  butalso  for  single-crystal  dielectric  plates  [3G]  -[38], 
Although  the  surface  wave  propagation  losses  in  amorphous  materials  become 


excessive  for  frequencies  above  100  MHz,  single-crystal  dielectric  plates  have 
been  used  at  frequencies  above  2 GHz.  Surface  waves  are  nondispersive  only  on 
the  surface  of  a semi-infinite  solid,  but  particle  motion  diminishes  exponentially 
with  distance  away  from  the  surface,  and  propagation  on  a plate  only  a few  wave- 
lengths thick  is  essentially  nondispersive.  Surface  waves  on  plates  more  than  a 
few  wavelengths  thick  are  the  slowest  of  the  nondispersive  waves,  having  veloci- 
ties between  0.87  and  0.96  of  the  shear  wave  velocity.  Thus,  for  a fixed  ITV 
product  and  \VDoX , an  implementation  of  a UTF  using  surface  waves  as  the  mode 
of  propagation  will  be  smaller  than  an  implementation  using  any  other  mode. 

Surface  waves  can  be  guided  in  several  ways  [39],  of  which  two  are  especially 
suitable  for  UTF  implementations : a thin  strip  on  the  surface  of  a plate,  with 
the  propagation  velocity  of  the  strip  lower  than  that  of  the  plate;  and  a thin  layer 
on  a plate  with  a slot  through  the  layer  to  the  plate,  with  the  propagation  velocity 
of  the  layer  greater  than  that  of  the  plate.  Both  of  these  guides  arc  dispersive. 
The  slot  guide  is  the  more  attractive  of  the  two,  because  it  has  the  least  disper- 
sion. There  is  a close  analogy  between  guided  surface  waves  and  guided  micro- 
wave  electromagnetic  waves  on  microstrips.  Microstrip  waveguide  components 
such  as  directional  couplers  and  hybrids  have  direct  surface  wave  equivalents  [-10], 
Surface  waves  have  velocity  and  wavelength  10s  times  smaller  than  electromag- 
netic waves  with  the  same  frequency.  Consequently,  surface  waveguide  compo- 
nents arc  approximately  105  times  smaller  than  their  microstrip  counterparts  or 
conversely,  for  a fixed  physical  size,  they  can  contain  10"'  times  more  compo- 
nents or  be  10?  times  more  complicated. 
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Amorphous  and  single-crystal  dielectric  materials  are  also  used  for  the 
propagation  of  bulk  waves,  both  longitudinal  and  shear,  for  frequencies  from 
about  1 MHz  to  100  MHz.  As  with  surface  waves,  bulk  wave  propagation  losses 
become  excessive  in  amorphous  materials  for  frequencies  much  greater  than 
100  MHz  and  single-crystal  materials  have  been  used  for  frequencies  in  excess 
of  2 GHz.  At  the  higher  frequencies,  the  wavelength  is  only  a few  microns  and 
well-collimated  acoustic  beams  with  diameters  on  the  order  of  1 millimeter  can 
be  generated  and  propagated  with  essentially  no  dispersion,  so  long  as  the  trans- 
verse dimensions  of  the  medium  are  sufficiently  larger  than  those  of  the  beam. 

In  general,  single-crystal  materials  are  elastically  anisotropic.  Therefore, 
the  orientation  of  the  crystal  axes  with  respect  to  the  direction  of  propagation  is 
an  important  consideration,  since  only  in  special  directions  will  the  direction  of 
energy  flow  be  colinear  with  the  normal  to  the  wavefront  [41]. 

Another  aspect  of  propagation  which  is  important  in  delay  line  implementa- 
tions is  the  propagation  delay  time  stability,  especially  with  respect  to  tempera- 
ture changes.  The  significance  of  delay  time  changes  for  the  tracking  of  a UTF 
matched  filter  was  discussed  earlier.  Delay  time  stability  is  also  important  in 
recirculating  time  compressors  such  as  the  DELTIC.  Ilorologic  steels  have 
been  developed  which,  in  addition  to  being  isopaustic,  are  magnetostrictivc  [42]. 
Wires  made  from  these  steels  propagate  torsional  mode  waves  with  good  delay 
time  stability  and  low  propagation  loss  and  can  be  tapped  magnetostrictively. 
Although  isopaustic  glasses  have  been  developed,  they  have  large  propagation 
losses  relative  to  fused  silica  [43];  however,  fused  silica  has  a large  temperature 


38 


coefficient  of  delay  and  must  be  temperature  controlled.  Single-crystal  quartz 


is  an  outstanding  material  in  this  context,  even  though  its  propagation  losses  at 
high  frequencies  are  somewhat  greater  than  for  other  crystalline  materials  in 
common  use,  since  delay  stable  orientations  of  the  crystal  axes  exist  for  shear 
waves  both  at  room  and  elevated  temperatures.  Because  of  the  importance  of 
time  delay  stability  to  microwave  acoustic  devices,  a considerable  expenditure 
of  effort  could  be  justified  in  the  search  for  additional  crystalline  material  with 
orientations  that  provide  delay  stability. 

Although  many  types  of  interactions  exist  between  acoustic  waves  and  other 
forms  of  energy,  the  most  commonly  used  types  for  the  transduction  of  energy 
between  electrical  and  acoustic  forms  are  the  electromechanical  and  magneto- 
mechanical  interactions.  Magnctostrictive,  ferroelectric  and  piezoelectric 
transducers  are  particularly  important  for  UTF  implementations.  These  trans- 
ducers are  reciprocal  and  can  be  used  both  to  launch  and  to  tap  a wave.  Another 
type  — the  optoacoustic  interaction  — is  gaining  importance  as  a means  of 
tapping  acoustic  waves.  In  general,  this  type  is  not  reciprocal  and  can  be  used 
only  for  tapping  the  acoustic  beam. 

Transduction  in  magnetos trictive  materials  is  similar  to  that  in  ceramic 
ferroelectric  materials.  In  both  types,  a bias  field  is  required  for  linear  opera- 
tion, and  a permanent  field  in  the  material  can  provide  this  bias.  Some  of  the 
magnctostrictive  materials  and  all  of  the  ferroelectric  ceramics  conmonly  used 
for  transduction  can  be  self-biased. 
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There  are  two  fundamental  types  of  tran  lion  in  these  materials: 


1 


dilational,  when  the  signal  field  is  colinear  with  the  bias  field;  and  rotational, 
when  the  signal  field  is  orthogonal  to  the  bias  field  [11].  The  colinear  trans- 
ducer is  used  for  longitudinal  and  extensions]  waves  and  the  crossed-field  trans- 
ducer is  used  for  shear  and  torsional  waves.  Both  magnetostrictive  and  ferro- 
electric materials  ai’e  used  extensively  for  these  transducers.  Both  types  of 
materials  are  used  in  colinear  transducers  for  extensions!  waves  in  wires  and 
ribbons  and  in  crossed-field  transducers  for  ■ ! ir  waves  in  strips  and  plates 
and  torsional  waves  in  wires.  When  the  mod.  , itself  is  magnetostrictive,  these 
various  modes  can  be  launched  and  tapped  di>  .-My  in  the  medium.  Ceramic 
transducers  can  be  bonded  directly  to  the  edi,  r of  ribbons,  strips  and  plates 
with  some  possible  advantages  relative  to  exit"  w!  magnetostrictive  transducers 
in  such  applications.  However,  a good  imped  ee  match,  wliich  is  sometimes 
difficult  to  obtain,  must  exist  between  the  bond.:.!  transducer  and  the  line;  other- 
wise the  transducer  may  have  small  fractional  bandwidth  [45]. 

Satisfactory  theories  for  the  operation  of  those  transducers  for  all  of  the 
propagation  modes  involving  linear  particle  displacement  wore  developed  many 
years  ago.  However,  until  a little  more  than  10  years  ago,  the  various  theories 
proposed  for  torsional  transducers  with  circular  cross-sections  were  in  serious 
disagreement  with  experimental  observation;.  [IS].  In  1958,  Yamamoto  [47]  pro- 
vided an  analytical  description  of  such  Ira  ns  d c.  •?,  which  is  in  excellent  agree- 
ment with  experimental  observations.  Allhou.  h the  analysis  was  applied  to 
magnetostrictive  transducers,  it  can  be  exit"  1-d  to  the  ferroelectric  case. 
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Single-crystal  piezoelectric  and  ferroelectric  transducers  arc  used  for  all 


modes  of  propagation  and  at  all  frequencies,  from  subaudio  through  microwaves. 
At  frequencies  above  100  MHz  the  control  of  crystal  thickness  mid  satisfactory 
bonding  techniques  become  critical  problems  [48].  The  theory  of  these  trans- 
ducers is  complicated  because  of  their  anisotropic  electric  and  elastic  proper- 
ties. The  search  for  new  materials  with  improved  properties  is  currently  an 
active  field. 

Thin-film  piezoelectric  transducers,  vacuum-deposited  on  the  propagation 
medium,  are  widely  used  for  longitudinal  and  shear  waves  at  frequencies  above 
500  MIIz  [49] -[51].  The  most  common  materials  are  cadmium  sulfide  and  zinc  ox- 
ide, with  the  latter  providing  somewhat  better  insertion  loss  and  bandwidth  charac- 
teristics than  the  former.  This  advantage  is  offset  to  some  extent  by  the  fact 
that  the  vacuum  deposition  of  zinc  oxide  is  more  difficult  than  that  of  cadmium 
sulfide.  For  both  materials,  the  best  performance  is  obtained  by  the  application 
of  thin-film  impedance  matching  layers  between  the  transducer  and  medium.  As 
with  the  single  crystal  piezo-  and  ferroelectric  materials,  new  materials  with 
improved  properties  for  thin-film  transducers  are  being  actively  sought. 

Optoacoustic  interactions  are  used  for  the  tapping  of  acoustic  waves  over  a 
wide  range  of  frequencies,  from  less  than  1 MHz  to  more  than  2 GHz.  The  tap- 
ping can  be  accomplished  by  the  modulation  of  a light  beam  incident  on  a wave 
which  is  interior  to  or  on  the  surface  of  an  optically  transparent  acoustic  medium. 
The  light  is  spatially  modulated  by  the  acoustic  wave,  according  to  the  phase  and 
amplitude  of  the  various  spectral  components  in  the  acoustic  wave  and  can  be 
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regarded  as  a continuously  distributed  tap  array.  Since  neither  the  bulk  wave 


nor  the  surface  wave  tap  is  generally  reciprocal,  they  are  not  used  for  launch 
transduction.  The  theory  of  diffraction  can  be  used  to  describe  the  behavior  of 
the  tap  array  and  the  resulting  analysis  usually  depends  on  the  detailed  nature 
of  the  optical  and  acoustic  arrangement. 

Magnetostrictive  Filters 

The  torsional  magnetostrictive  wire  delay  line  is  unusually  simple  in  con- 
cept and  implementation.  One  or  more  turns  of  wire  encircling  the  delay  line 
can  serve  as  either  a launch  transducer  or  a tap  with  the  weight  controlled  by 
the  number  of  turns  and  their  direction.  The  magnetic  bias  required  for  linear 
magnetostrictive  transduction  is  obtained  by  passing  a current  down  the  delay 
line  wire,  thus  establishing  a circumferential  crosscd-ficld  bias  within  the  de- 
lay line.  If  the  delay  line  material  is  sufficiently  remanent,  the  current  can  be 
removed  and  a permanent  bias  field  will  remain.  When  a current  is  passed 
through  a launch  wire,  a longitudinal  magnetic  field  is  established  in  the  delay 
line  which  acts  as  the  signal  field.  The  bias  field  and  signal  field  interact  with 
the  magnetostrictive  properties  of  the  delay  line  to  create  a torsional  strain 
which  propagates  along  the  line  (the  Wiedemann  effect  [52]).  Conversely,  a tor- 
sional strain  in  the  delay  line  interacts  with  the  bias  field  and  the  magnetostric- 
tion to  create  an  additional  longitudinal  field  component  (the  inverse  Wiedemann 
effect)  and  as  the  strain  propagates  along  the  lino  and  passes  a tap,  this  asso- 
ciated field  induces  an  cmf  in  the  tap.  Yamamoto  gives  the  relationship  between 
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the  torsion  per  unit  length,  0,  the  longitudinal  mager-fostrictivc  coefficient,  X,, 
and  the  magnetic  field  intensities,  IIr,  H,.,  and  11,,  as 

0 - (19) 

where  11®  - H®  + Il£,  IIX  and  IIC  are  the  longitudinal  and  circumferential  field  in- 
tensities, respectively,  and  r is  the  radius  of  the  d lay  line  wire  [47].  (19)  indi- 

cates that  when  H:  « IIC,  the  torsion  is  linear  in  1,. 

Torsional  delay  line  UTFs  have  been  dovcloi  ’ and  fabricated  by  Whilehouse 
and  Lindsay,  operating  at  center  frequencies  in  !'  region  of  1 MHz,  with  almost 
100%  fractional  bandwidth  and  with  T\V  products  c ! pproximatcly  10"  [25].  The 
T\V  product  of  these  filters  can  be  increased  to  approximately  101  by  compensat- 
ing the  propagation  losses,  and  it  may  be  possibl  ! . increase  the  operating  fre- 
quency  to  about  10  MHz  without  decreasing  the  f ■ '<  a:  1 bandwidth,  by  employ- 
ing a delay  line  wire  of  smaller  diameter  and  rod- • h ring  the  tap  structure.  How- 
ever, advances  beyond  these  values  are  blocked  by  the  fundamental  limitations 
due  to  the  polycrystalline  nature  of  the  delay  line  wire. 

The  torsional  delay  line  UTF  is  especially  : i!  dde  for  the  implementation  of 
multiple -port  arrays.  In  one  implementation,  a single  turn  tap  is  passed  through 
a small  linear  ferrite  core  and  serves  not  only  a , r tap,  but  also  as  the  primary 
winding  of  a wide-band  transformer.  The  weigl  ( of  I he  tap  is  determined  by  the 
direction  and  number  of  the  secondary  turns  on  th  core  and  the  summation  is 
accomplished  by  connecting  the  secondary  windii  g ; of  all  the  taps  in  series.  If 
the  core  at  each  tap  has  a large  enough  aperture  In  p i mil  a number  of  secondary 

43 


windings,  then  a corresponding  numb e ••  of  independent  weighting  functions  and 
multiple  ports  can  be  obtained,  hi  actual  practice,  such  a filter  might  be  con- 


structed by  passing  a shuttle,  carrying  a continuous  wire,  through  each  core 
in  turn,  in  a manner  analogous  to  the  si',  lung  operation  in  sewing.  The  tap 
weights  are  determined  by  the  direct!  a d number  of  times  the  shuttle  is 
stitched  through  each  core  before  pass!  >g  on  to  the  next  one.  Each  time  the 
shuttle  is  stitched  through  the  complet  sot  of  cores,  it  carries  a wire,  the  two 
ends  of  which  constitute  a port.  Thus.,  Ih  > number  of  independent  ports  possible 
is  limited  by  the  number  of  wires  11 1;  . < .,  > be  passed  through  the  core  apertures 
without  blocking  the  passage  of  the  shnbV.-  A single  port  filter  of  this  type  is 
shown  schematically  in  Fig.  9 and  a section  of  an  eight-port  filter  operating  at 
a center  frequency  of  1 MHz,  with  la  r fractional  bandwidth  capability,  and 
with  a maximum  TW  product  of  512  i ’ own  in  Fig.  10  [25]. 

The  number  of  possible  indepenu  1 ' ports  can  be  increased  by  using  split 
cores.  In  this  case,  the  secondary  wi>  Vugs  arc  prefabricated  from  a continuous 
wire  and  then  slipped  over  one  scctii  . of  each  of  the  split  cores,  after  which  the 
other  core  sections  arc  placed  in  posith  < to  close  the  cores.  Such  a multiple 
port  filter  is  shown  schematically  in  1 ■.  11.  The  fabrication  process  for  the 
secondary  windings  can  be  performed  or.  a loom.  The  warp  would  consist  of  2n 
wires,  where  n is  the  maximum  weight  for  any  tap  in  the  filter,  half  of  the  wires 
for  positive  weights  and  half  for  negati  ve  weights.  At  each  position  along  the 
wires  corresponding  to  the  tap  positio.  , |w;|  wires  arc  lifted  — |\v,|  from  the 
positive  group  for  a positive  weight,  < v [\Vjj  from  the  negative  group  for  a 
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negative  weight,  where  w,  is  the  weight  of  the  i tap.  At  each  tap  position, 
then,  a number  of  wires  corresponding  to  the  magnitude  of  the  weigh!  at  that  tap 
position  has  been  lifted  from  the  group  of  wires  corresponding  to  the  sign  of  the 
weight.  The  warp  is  then  gathered  into  a compact  bundle  of  wires,  care  being 
exercised  to  keep  the  appropriate  wires  lifted.  The  split  core  halves  are  then 
inserted  into  the  spaces  between  the  lifted  and  undisturbed  wires  and  11  other 
halves  are  used  to  close  the  cores.  Finally  all  of  the  wires  are  com  d in 
series,  with  the  negative  wires  bucking  the  positive  wires.  A tapping  structure 
results,  with  the  correct  weights.  For  multiple  ports,  2nm  wires  are  used, 
where  m is  the  number  of  ports,  with  the  wires  connected  in  groups  c ! 2n  for 
each  port.  A Jacquard  loom,  such  as  the  one  used  at  the  Massnehus.  ' s 
InstiUitc  of  Technology  for  the  fabrication  of  "braid"  memories,  can  1 used  to 
loom  these  tapping  structures  [53]. 

The  torsional  delay  line  UTF  described  above  has  fixed  weighting;  f actions. 
Another  type  of  wire  delay  line  UTF,  using  extcnsional  waves,  has  an  a lterable 
weighting  function.  This  filter  makes  use  of  a nondestructive-read,  ferro- 
electric memory  [54]  which  uses  the  coincidence  between  a high  intensity  exten- 
sional  wave  pulse  propagating  along  the  delay  line  and  a time  varying  current  in 
the  delay  line  wire  to  form  a remanent  field  in  the  delay  line  material.  The  mag- 
nitude and  sign  of  the  remanent  field  at  a particular  point  on  the  line  ai>  con- 
trolled by  the  magnitude  and  sign  of  the  current  in  the  delay  line  wire  at  the  time 
the  acoustic  pulse  passes  the  point.  An  acoustic  pulse  launched  along  th  line 
induces  an  cmf  in  the  delay  line  wire  given  by  h(t)  w(t),  where  v (z/c)  i the 
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function  describing  the  remanent  field.  Thus,  the  voltage  appearing  across  the 
ends  of  the  delay  line  wire  is  the  weighted  and  summed  output  from  contiguous 
taps,  with  the  weighting  function  w(t).  To  replace  the  weighting  function  with  a 
new  one,  a time  varying  current  of  the  desired  form  is  passed  through  the  wire, 
and  an  acoustic  pulse  is  launched  along  the  delay  line.  A UTF  of  this  type  can 
be  made  to  operate  at  a frequency  of  1 MHz,  and  100%  bandwidth  with  a TW  prod- 
uct of  lO3  might  be  possible.  The  TW  product  in  this  filter  is  limited  by  the  fact 
that  the  high  intensity  acoustic  pulse  needed  to  set  the  remanent  field  drives  the 
material  of  the  delay  line  into  a nonlinear  acoustic  region,  with  large  propagation 
losses  as  a consequence. 

Another  ferro-acoustic  delay  line,  originally  developed  as  a stress-current 
coincidence  memory  [55],  has  interesting  possibilities  as  a delay  line  for  UTFs 
when  used  in  a current-current  coincidence  mode  similar  to  that  used  for  plated 
wire  memories.  The  delay  line  consists  of  a beryllium -copper  wire,  plated 
with  a thin  permalloy  film.  The  remanent  filed  is  set  by  the  coincidence  of  a 
current  in  the  delay  line  wire  and  a current  in  a wire  perpendicular  to  the  delay 
line.  The  impulse  response  is  related  to  the  rcmanence  by  the  same  equation  as 
for  the  stress-current  coincidence  filter:  h(t)  w(t).  Since  in  this  filler  the 
rcmanence  is  not  set  acoustically,  the  limitations  on  the  TW  product  caused  by 
the  rapid  attenuation  of  the  acoustic  setting  pulse  are  absent  and  it  should  be 
possible  to  construct  filters  operating  at  10  MHz,  with  100%  bandwidth  and  with 
TW  products  of  approximately  104. 
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A UTF  with  multiple  output  ports  as  well  as  alterable  weighting  functions  is 


possible  if  the  wire  delay  line  is  replaced  by  a glass  or  fused  silica  plate  with  a 
vacuum-deposited  thin  magnetostrictive  film  on  its  surface.  As  with  the  wire 
delay  lines,  the  state  of  remanence  can  be  set  by  either  stress-current  or 
current-current  coincidence.  If  the  acoustic  transducer  is  many  wavelengths 
across,  the  acoustic  wave  will  be  well  collimated  and  plane  wave  propagation 
will  occur.  Thus,  the  acoustic  stress  will  be  constant  in  each  differential  strip 
transverse  to  the  direction  of  propagation.  If  a number  of  wires  parallel  to  the 
direction  of  propagation  are  placed  in  juxtaposition  with  the  surface,  each  of 
these  wires  can  be  used  as  an  independent  channel  for  setting  the  remanence  in 
that  part  of  the  film  directly  under  the  wire.  Once  the  remanence  has  been  set, 
each  wire  becomes  an  independent  output  port.  The  number  of  ports  possible  is 
limited  only  by  the  spacing  required  between  channels  to  prevent  cross-talk,  and 
by  the  width  of  the  plate.  A similar  filter,  in  which  the  magnetostrictive  film 
has  uniaxial  anisotropy,  limits  the  weighting  function  to  binary  values  [56]. 
Filters  operating  at  frequencies  as  high  as  100  Mil/.,  with  100%  fractional  band- 
width and  with  a TW  product  of  104  might  be  possible. 

Surface  Wave  Filters 

One  of  the  most  interesting  advances  in  microwave  acoustics  is  the  use  of 
surface  waves  on  piezoelectric  and  ferroelectric  substrates  [57].  Although  sur- 
face waves  can  be  used  at  frequencies  from  below  1 MHz  to  above  1 GIIz,  the 
most  promising  applications  for  UTFs  are  at  frequencies  from  about  10  MHz  to 
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500  MHz.  For  these  frequencies,  surface  wave  devices  are  small  in  size  and 


are  easily  fabricated  by  techniques  similar  to  those  used  by  the  semiconductor 
industry  for  large  scale  integrated  circuits. 

Many  of  the  piezoelectric  materials  used  as  substrates  are  also  semicon- 
ductors, and  amplification  of  surface  waves  on  such  materials  has  been  achieved 
by  the  interaction  of  the  electric  field  of  the  surface  wave  with  charge  carriers 
drifting  at  the  surface  wave  velocity  [58].  Since  these  fields  extend  above  the 
surface,  a nonsemiconducting  substrate  can  be  used  in  conjimction  with  a non- 
piezoelectric  semiconductor:  the  two  materials  need  only  to  be  brought  into  prox- 
imity with  a drift  field  applied  to  the  semiconductor  [59],  This  separation  of  the 
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functions  of  acoustic  propagation  and  semiconductor  amplification  allows  each 
material  to  be  optimized. 

Surface  waves  on  isotropic  nonpiezoelectric  substrates  have  been  studied 
extensively  [35].  Although  transduction  on  such  substrates  is  difficult,  there  is 
still  considerable  interest  in  them  due  to  their  superiority  for  guided  surface 
waves  [60].  The  most  successful  transducer  for  nonpie zoelecti'ic  substrates  is 
the  Sokolinskii  comb  transducer  [61].  This  transducer  consists  of  a longitudinal 
mode  piezoelectric  crystal,  bonded  to  one  face  of  a buffer  plate  whose  opposite 
face  has  periodic  grooves  and  is  in  contact  with  the  substrate.  Electrical  exci- 
tation of  the  crystal  produces  a periodic  stress  distribution  on  the  surface  in  con- 
tact with  the  comb,  due  to  the  grooves  in  the  plate.  This  technique  results  in  an 
efficient  surface  wave  transducer  with  good  directivity,  and  with  the  direction  of 
surface  wave  propagation  normal  to  the  grooves. 


. 
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For  piezoelectric  materials  the  Wiiite-Voltmcr  interdigitated  electrode 


transducer  has  been  highly  successful.  This  transducer  consists  of  an  array  of 
periodically-spaced  parallel  conductive  strips  on  the  surface  of  the  substrate, 
with  alternate  strips  electrically  connected  together  to  form  a comb-like  array 
of  interdigitated  electrodes  [62].  When  the  array  is  excited,  periodic  stx-esses 
are  generated  at  the  surface  by  the  interaction  of  the  resulting  electric  field 

1 

with  the  piezoelectric  material,  and  a surface  wave  is  launched. 

The  Sokolinskii  transducer  must  have  a large  number  of  grooves  for  low  in- 
sertion loss  and  has  an  imp  1 ■ response  which  is  approximately  a finite  duration 
sine  wave,  with  the  number  of  cycles  equal  to  the  number  of  grooves.  Conse- 
quently,  it  is  a narrowband  transducer.  The  White- Voltmer  transducer  requires 
a large  number  of  strips  for  low  insertion  loss,  when  the  piezoelectric  coupling 
coefficient  is  small,  and  h . an  impulse  response  similar  to  that  for  the 
Sokilinskii  transducer.  Thus  it  is  also  a narrow  band  transducer.  However, 
when  the  coupling  coefficient  is  large,  low  insertion  loss  can  be  achieved  with  a 
few  electrodes,  and  in  this  case  the  transduction  is  wide  band  [63]. 

The  surface  wave  tram  lucer  array  is  an  example  of  a distributed  array  of 
the  type  described  in  the  discussion  of  generalized  transversal  filters.  If  the 
input  array  has  weighting  fi  a t ion  wt  and  the  output  array  has  weighting  function 
w0,  then  the  impulse  response  of  the  array  pair  is  wf  *w0  which  is  the  cross- 
correlation of  the  input  and  (!.?  output  array  weighting  functions.  Transducer 
arrays  of  the  types  describ  d above,  when  used  as  input  and  output  arrays,  have 
an  impulse  response  (the  ere:  s correlation  of  two  sine  waves  of  equal  duration) 
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which  is  a sinusoid  with  a triangular  envelope  with  a duration  twice  that  of  the 


impulse  response  of  the  individual  arrays.  Thus  an  array  pair  acts  as  a narrow 
band  filter  whose  bandwidth,  for  a fixed  array  length,  is  inversely  proportional 
to  the  number  of  strips  or  grooves  and  whose  frequenc}'  transfer  function  is  of 
the  form  sin?x/xs. 

The  result  indicates  that  spectrum  analyzer  can  be  constructed  with  high 
resolution  and  low  si  delobe  response.  If  more  than  one  pair  of  such  electrode 
arrays  are  cascaded,  then  the  total  frequency  transfer  function  will  have  the  form 
sin8rx/xSr,  where  r is  the  number  of  pairs  in  cascade.  For  r = l the  first  minor 
lobe  is  down  2G  dn,  and  for  r 2 the  first  minor  lobe  is  down  more  than  50  dB. 

In  addition,  when  two  sets  of  such  cascaded  transducer  pairs  are  deposited  on 
the  same  crystal  substrate  and  driven  from  separate  sources,  the  phase  differ- 
ence between  the  two  signals  can  be  measured  simultaneously  with  their  spectra 
[64]  since  differential  phase  stability  is  assured  by  the  use  of  a common  substrate. 

In  many  applications  it  is  desirable  to  have  frequency  transfer  functions 
w'hich  are  either  broad  band  or  have  small  side  lobes.  Two  types  of  weighting 
functions  which  might  be  used  to  achieve  these  responses  are  "period  variation" 
weighting,  where  the  spacing  between  interdigitations  changes  along  the  array, 
and  equally-spaced  interdigitations  with  amplitude  modulation  of  the  individual 
tap  elements.  An  essentially  rectangular  wide  band  frequency  response  can  be 
achieved  with  the  former,  if  the  impulse  response  of  the  array  has  an  instantan- 
eous frequency  which  changes  linearly  with  time;  a sidelobe-free  frequency  re- 
sponse can  be  achieved  with  the  latter,  if  the  amplitude  modulation  is  Gaussian. 
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A third  type  of  weighting  function,  in  which  the  intei’digitated  transducer 


elements  arc  phase  modulated,  is  required  for  many  applications.  An  example 
of  a phase-modulated  code  with  a narrow  correlation  function  and  with  uniformly 
small  sidelobes  is  a Barker  code  [65].  A transducer  array  using  simple  elcc- 
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trodc  pair  inversion  to  implement  a 7 bit  Barker  code  is  shown  in  Fig.  12a, 

f 

along  with  the  corresponding  electric  fields.  Electrode  pair  inversion,  however, 

1 

fails  to  provide  the  impulse  response  expected  of  a Barker  code,  because  of  the 
cquipotential  regions  formed  at  the  phase  changes  of  the  code,  with  consequent 
loss  of  spatial  components  in  the  stress  distribution. 

A phase-coded  transducer  which  has  the  required  stress  distribution  has 
been  developed  by  Whit  chouse  and  Lindsay,  and  utilizes  one  additional  electrode 
with  a reconnection  of  the  other  electrodes.7  This  transducer  and  its  electric 
fields  are-shown  for  comparison  in  Fig.  12b.  By  appropriately  weighting  the 
amplitude  of  a transducer  composed  of  these  phase -invertible  elements,  an  arbi- 
trary impulse  response  can  be  synthesized  by  the  techniques  given  in  the  first 
section  of  this  paper. 

As  an  example  of  the  synthesis  of  an  arbitrary  impulse  response,  a pulse-in 
to  pulse-out  distributed  surface-wave  transducer  delay  line  is  described.  It  is 
possible  to  obtain  pulse-in  to  pulse-out  behavior  if  the  input  and  output  arrays  are 
each  a simple  electrode  pair.  However,  a large  number  of  elements  in  each 
array  is  required  to  achieve  good  signal-to-noise  ratios  on  substrates  with  low 
electromechanical  coupling  coefficients.  Arrays  with  large  numbers  of  elements 

7 

A paper  describing  tin’s  work  in  more  detail  is  in  preparation. 
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can  be  designed  for  pulse-in  to  pulse-out  operation,  but  care  must  be  taken  to 


choose  array  weightings  whose  corrlcation  functions  have  low  sidelobes  in  order 

to  avoid  signal  interference  when  continuous  signals  are  transmitted. 

« 

The  Huffmann  codes  [GG]  are  phase-reversal  codes  with  impulse-equivalent 
behavior.  Thcjr  exist  for  all  code  lengths,  and  have  only  one  small  correlation 
sidelobe  on  each  side  of  the  central  impulse.  However  these  codes  are  not 
attractive  for  surface  wave  delay  lines  because  they  require  a large  variation  of 
the  array  element  amplitudes,  with  a consequent  difficulty  of  fabrication.  "I'll i s 
problem  can  be  alleviated  by  the  use  of  a pair  of  binary  phase-reversal  c odes 

' 

which  have  the  property  that  the  sum  of  the  correlation  functions  of  the  code 

• < 

pairs  is  an  impulse,  with  no  side  lobes.  Golay  has  studied  such  codes,  which  lie 

% 

calls  "co mp  1 i mentnry  scries.  " For  their  properties,  the  reader  is  referred  to 
his  work  [67].  Since  these  codes  arc  phase  reversal  codes,  the  three-electrode 
differential  transducer  is  a natural  choice.  In  phase-coded  arrays  that  have 
these  complimentary  scries  as  weighting  functions,  the  number  of  positive  elec- 
trodes is  not  equal  to  the  number  of  negative  electrodes.  Consequently,  the 
transducer  is  capacitivcly  imbalanced.  However,  the  correlation  function  of  the 
negative  of  a function  is  the  same  as  that  of  the  function  itself.  Thus,  one  of  the 
pair  of  codes  can  be  negated,  so  that  a parallel  connection  of  the  two  transducers 
will  have  capacitive  balance.  When  the  weighting  functions  arc  so  chosen,  the 
transducer -can  be  conveniently  driven  by  a wideband  differential  transformer, 
with  some  gain  in  rejection  of  interference.  A delay  line  using  a 128-element 
Golay  complementary  series  array  connected  for  capacitive  balance  is  shown 
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in  Fig.  13.  The  implementation  of  complementary  series  phase-coded  surface- 
wave  transducers  should  mate  ilde  wide  band  pulse  delay  lines,  with  negligi- 
ble side  lobe  interference,  for  o;  ■ 'lion  to  300  Mil/..  'Ilicsc  delay  lines  should 
have  TW  products  of  about  10  ' a <1  they  can  be  fabricated  with  existing  photo- 
resist technology.  These  dcln\  lb  when  combined  either  with  microwave 
acoustic  amplification  or  mi  cm  •electronic  amplification,  can  be  used  for  trans- 
versal filter  time  compressors  or  for  DELTiCs,  as  well  as  for  matched  filters 
and  spectrum  analyzers. 


Optoacoustic  Filters 

With  the  recent  advent  of  lam  .*« , interest  in  the  optoacoustic  interaction  as 
a means  of  tapping  the  inform;  'in,  contained  in  an  acoustic  delay  line  has  in- 
creased. Typically,  the  inforr  : Lion  in  the  acoustic  beam  will  be  contained  in  a 
band  of  frequencies  about  a ce  1 carrier  frequency,  f0.  From  spatial  samp- 
ling considerations,  the  acoustic  beam  can  be  considered  to  be  composed  of 
contiguous  segments  of  length  a.,/2  c/2fa,  where  e is  the  acoustic  propagation 

velocity,  fB  is  the  largest  frequency  by  which  the  carrier  is  modulated,  and 

. 

fB  £ f0.  Each  of  these  segment:  will  interact  with  that  portion  of  the  incident 
light  beam  which  it  intercepts,  if  (lie  incident  light  beam  intercepts  a length 
A^/2  of  the  acoustic  beam,  then  tlv-ro  is  only  one  tap  on  the  delay  line,  and  the 
output  of  that  tap  will  simply  be  modulated  in  time  as  the  acoustic  waves  propa- 
gate by.  If  the  incident  light  bourn  intercepts  a length  larger  than  Ae/2,  then 
each  of  several  contiguous  sc;,.  of  the  acoustic  beam  will  interact  separately 
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128  ELEMENT  PHASE  REVERSAL  ARRAYS 


Fig.  13.  Surface-wave  ultrasonic  transversal  filter  with  phase-invertible 
elements,  employing  Golay  complimentary  series  codes.  Upper 
and  lower  arrays  are  complimentary. 


with  that  portion  of  the  incident  light  beam  which  it  intercepts,  so  that  several 
contiguous  optoacoustic  taps  result,  each  of  which  can  be  separately  weighted  in 
amplitude,  phase,  or  polarization. 

If  the  acoustic  beam  is  a bulk  wave  propagating  in  a liquid  or  a solid,  the 
delay  experienced  by  any  given  light  ray  traveling  through  the  acoustic  beam  will 
depend  on  the  degree  to  which  the  density  along  its  path  has  been  increased  or 
decreased  by  the  alternating  compressions  and  rarefactions  of  the  acoustic  beam. 

Not  only  will  the  delay  along  the  optical  path  be  alternately  increased  and  decreased 
by  such  action,  but  also  will  the  ray  be  refracted  or  polarized  by  the  acoustic 
medium,  resulting  in  both  phase  and  amplitude  corrugations  of  the  optical  wave- 
front.  Such  corrugations  will  generate  spatial  interference  patterns  similar  to 
those  caused  by  a conventional  diffraction  grating. 

A grating-like  optical  interference  pattern  can  also  be  generated  using  acous- 
tic surface  waves.  Such  waves  can  change  optical  path  delay  in  a manner  analogous 
to  the  bulk  wave  interaction  just  described  [GS ] , and  they  can  also  be  used  to  frus- 
trate the  total  internal  reflection  of  a light  beam  at  a solid-air  interface.  For  the 
latter  effect,  a second  solid  is  placed  adjacent  to  the  first,  and  the  surface  wave  is 
propagated  on  cither  of  the  two  adjacent  surfaces  so  that  a maximum  frustration 
occurs  at  the  crests  of  the  wave,  and  a minimum  frustration  occurs  at  its  troughs. 

This  gives  rise  to  a reflected  and  a transmitted  light  beam  both  of  which  arc 
modulated  in  amplitude  across  their  widths. 

In  many  cases  of  interest  only  the  light  diffracted  into  the  first  order  lobc(s) 

is  important,  and  the  light  in  other  orders,  including  the  zeroth  order,  can  be 
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neglected  or  can  be  eliminated  by  means  of  spatial  filters  or  polarization  analy- 


zers. In  general,  the  rays  diffracted  into  the  first  order  lobe  will  travel  at  some 
angle  relative  to  the  acoustic  wavefronts  and  thus  will  have  the  possibility  of  over- 
lapping with  rays  being  diffracted  by  an  optical  tap  contiguous  (o  (he  one  in  which 
the  ray  originated.  The  proportion  of  overlapping  rays  decreases  as  the  width,  d, 
of  the  acoustic  beam  is  decreased,  or  as  the  maximum  modulation  frequency  fB 
is  decreased.  For  acoustic  bulk  waves,  decreasing  the  acoustic  beam  width  will 
usually  cause  a decrease  in  the  amount  of  light  diffracted,  and  decreasing  fB 
reduces  the  bandwidth,  so  that  a tradeoff  between  acoustic  beam  width,  band- 
width, and  signal-to-noise  ratio  may  be  required.  Large  fractional  bandwidths 
on  the  order  of  50-100%  are  not,  in  fact,  possible  when  d » A^/X,  where  A0  ; c/f0, 
and  X is  the  optical  wavelength  (this  is  the  region  where  Bragg  diffraction  domi- 
nates Raman-Nath  diffraction  [GO] ).  The  interaction  of  light  with  surface  waves 
is  inherently  broadband,  since  for  surface  waves,  the  transverse  dimension  d is 
of  the  order  of  A0. 

An  optical  system  can  be  regarded  as  a two-dimensional  analog  of  an  elec- 
trical system,  with  the  input  and  output  planes  of  the  optical  system  correspond- 
ing to  the  input  and  output  ports  of  the  electrical  system.  It  is  known  that  the 
optical  counterparts  of  various  electrical  filtering  operations  (c.g.,  low  pass, 
high  pass,  bandpass,  band  stop,  differentiation,  integration)  can  be  performed 
by  means  of  lenses  and  apertures  [70].  The  information  to  be  filtered  is  presented 
at  the  input  plane  as  a transmission  or  reflection  function  with  an  arbitrary  distri- 
bution of  amplitude,  phase  and  polarization,  and  the  input  plane  is  illuminated 
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with  coherent,  monochromatic  light.  Spatial  multiplication  occurs  when  a second 
transmission  or  reflection  function  is  placed  either  in  the  same  plane  as  that 
occupied  by  the  first  function,  or  at  an  image  of  that  plane.  An  optical  multiplier 
based  on  these  principles  is  illustrated  in  Fig.  14. 

When  the  transmission  function  moves  across  the  input  plane,  then  the  result- 
ing distribution  in  the  output  plane  will  correspondingly  change  in  time  in  a man- 
ner depending  upon  the  intervening  optical  system.  If  the  intervening  optical 
system  serves  to  spatially  multiply  the  input  transmission  function  by  a second 
transmission  function  (the  weighting  function),  then  the  resulting  spatial  product 
is  equivalent  to  a weighted  tapping  of  the  input  transmission  function.  Optical 
integration  of  this  weighted  tapping  forms  an  optical  transversal  filter. 

If  the  optical  integration  is  astigmatic  so  (hat  it  integrates  this  product  only 
along  the  direction  of  motion,  then  a multiplicity  of  simultaneous  transversal 
filters  can  be  formed,  hi  particular,  if  the  input  plane  is  divided  into  a number 
of  strips  parallel  to  the  direction  of  motion,  then  each  strip  and  its  corresponding 
output  can  be  independent  of  all  other  strips  and  outputs. 

An  example  of  an  optical  transversal  filter  is  a film  strip  moving  along  an 
axis  in  the  input  plane.  Unfortunately,  film  processing  delays  and  slow  film 
movement  make  this  form  of  optical  filter  too  slow  for  high-speed  real-time 
signal  pi'oeessing.  High-speed  and  real-time  capability,  however,  are  available 
in  ultrasonic  delay  lines  in  which  light  interacts  photoelastically  with  the  acoustic 
waves  in  an  optically-transparcnt  medium.  Although  ultrasonic  delay  lines  do  not 
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have  the  high  two-dimensional  information  storage  capability  of  film,  a large 
number  of  them  can  be  used  side-by-side,  filling  up  the  whole  input  plane,  to 


partially  offset  this  loss. 

Figure  15  shows  two  examples  of  such  an  optoacoustic  filter  in  which  the 
electrical  signal  of  interest  has  been  transformed  into  an  acoustic  signal  to  be 
interrogated  by  the  light  beam  from  the  left.  The  modulated  light  beam  is  then 
weighted  by  (or  compared  to)  a reference  function  as  represented  by  a transpar- 
ency or  phase  plate  or  coded  polarization  sheet  (Fig.  15a),  or  by  another  ultra- 
sonic signal  in  a second  delay  line  (Fig.  15b),  and  finally  integrated  and  presented 
at  the  output  of  a light  detector.  This  second  representation  illustrates  the  ability 
of  the  optoacoustic  filter  to  program  the  reference  function  in  real  time.  The 
positions  of  the  signal  and  reference  planes  can  be  interchanged,  but  their  separa- 
tion (or  the  separation  between  one  plane  and  the  image  of  the  other)  must  be  small 
enough  to  prevent  diffraction  overlap  of  adjacent  optical  information  cells. 

In  Fig.  16  is  the  diagram  of  a device  utilizing  two  acoustic  cells  oriented  to 
take  advantage  of  Bragg-angle  incidence  of  the  optical  beams.  One  cell  is  placed 
in  juxtaposition  to  the  other  cell  (or  its  image)  to  re-diffract  the  first-order  lobe 
transmitted  through  the  polarizer  P,.  When  the  two  acoustic  signals  are  matched, 
the  net  phase  change  of  any  light  ray  undergoing  the  double  diffraction  is  constant, 
regardless  of  the  acoustic  frequency,  so  that  the  wave  passing  through  P2  is  plane. 
When  the  signals  are  not  matched,  the  optical  wavefront  passing  through  P2  may 
be  a plane  wave  of  lesser  intensity,  a plane  wave  in  a different  direction,  or  a 
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non-plane  wave,  so  that  the  intensity  at  the  focal  point  of  the  lens  is  below  the 
maximum  given  by  matched  signals.  A slightly  Doppler-shifted  version  of  the 
reference  signal  will  produce  an  output  plane  wave  which  will  focus  to  a different 
point  in  the  focal  plane  < ' the  lens.  If  the  polarizers  Px  and  P2  ai’e  deleted,  some 
other  means  (c.g.,  spat  i : 1 filtering,  or  optical  heterodyning)  must  be  used  to 
separate  the  twice  difir  1 beam  from  the  zcro-order(undiffracted)  light 
collinear  to  it.  Squire  mil  Alsup  [18]  have  constructed  a device  of  this  type 
which  can  operate  at  a fr  yeney  of  1 GHz  with  50%  bandwidth  and  a TW  product 
of  103. 

. * 

A modification  of  the  two-cell  Bragg  filter  described  above  is  presently  under 
investigation  and  provide  a simple  means  of  separating  the  undiffracted  light  from 
the  signal-carrying  diflh : < ' ; d light.  The  two  configurations  being  investigated  are 
shown  in  Fig.  17.  In  (h  first,  the  two  cells  are  placed  so  that  the  axes  of  acoustic 
pi'opagation  are  orient  i!  ; (he  angle  0-0,  where  0 is  the  Bragg  angle  for  a 
carrier  frequency  fx,  and  ^ is  the  Bragg  angle  for  a carrier  frequency  f2.  In  either 
configuration,  the  beam  f light  which  is  diffracted  by  an  unmodulated  carrier  fx 
in  the  first  cell  will  be  incident  at  (he  appropriate  angle  to  be  rcdiffracted  by  an 
unmodulated  carrier  f.  in  the  second  cell- in  the  direction  specified  by  the  Bragg 
conditions  for  the  two  ( .11s.  Furthermore,  an  unmodulated  carrier  frequency 
fx  + Af  in  the  first  cell  w ill  diffract  a Bragg  beam  which  will  be  rediffracted  by 
an  unmodulated  carrier  frequency  f2  H A f in  the  second  cell  into  the  same  direc- 
tion indicated  above  for  frequencies  fx  and  fc;  the  plus  sign  is  appropriate  to 
Fig.  17a,  and  the  minu.  to  Fig.  17b.  In  both  cases,  the  twice-diffracted  light 
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emerges  in  a direction  not  collinear  with  the  incident  light  (except  when  0 is  equal 


to  0 in  the  first  configuration),  so  that  the  undiffractcd  light  can  be  removed  by 
spatial  filtering  in  the  focal  plane  of  the  integrating  lens. 

If  one  of  the  acoustic  cells  in  the  preceding  example  is  replaced  by  a fixed 
phase  or  amplitude  grating  designed  to  represent  the  signal  to  be  matched,  then 
a device  results  such  as  that  shown  in  Fig.  18  [71],  where  the  loss  of  versatility 
provided  by  an  alterable  weighting  function  has  been  offset  by  a gain  in  light  signal 
intensity. 

A single  acoustic  cell  can  act  as  a Fresnel  lens  while  maintaining  a unique 
Bragg  diffraction  geometry  when  the  signal  to  be  processed  is  a linear  FM  sweep, 
or  ’’chii'p"  [72],  [73].  As  diagrammed  in  Fig.  19,  this  configuration  satisfies  the 
Bragg  diffraction  condition  simultaneously  for  all  frequencies  in  the  acoustic 
beam  at  only  one  instant  of  time,  and  it  is  at  that  instant  that  the  acoustic  beam 
focuses  the  diffracted  portion  of  the  diverging  input  light  to  a single  point  without 
the  use  of  auxiliary  lenses.  An  alternative  form  of  this  technique  uses  a col- 
limated incident  light  beam  and  anisotropic  polarization  rotation  of  the  focused 
diffracted  beam  [74].  It  should  be  possible  to  achieve  bandwidths  and  TW  products 
in  this  type  of  filter  comparable  to  those  given  in  the  two  previous  examples. 

The  frustrated  total  internal  reflection  by  means  of  surface  waves,  described 
earlier,  is  utilized  in  the  filter  shown  in  Fig.  20  [75].  The  beam  diffracted  from 
the  first  sui'face  wave  is  spatially-filtered  within  an  afocal  lens  system  which 
images  the  first  surface  wave  upon  the  second  surface  wave.  The  twice -modulated 
light  is  gathered  by  a lens  to  a point  at  which  is  placed  a pliotodctector,  whose 
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Optoacoustic  transversal  filter  utilizing  frustrated  total  internal 
reflection  of  light  by  surface  waves. 
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output,  represents  the  correlation,  or  convolution  of  the  two  acoustic  signals 
propagating  as  surface  waves.  Like  the  first  Bragg  device  described,  this  filter 
has  the  capability  of  an  alterable  weighting  function,  but  its  bandwidth  limitations 
are  governed  by  realizable  carrier  frequencies,  instead  of  beam  width.  Such  a 
filter  should  have  performance  characteristics  comparable  to  those  of  the  surface- 
wave  filters  described  earlier. 

CONCLUSIONS 

Linear  transversal  filters  arc  ideal  for  the  implementation  of  filters  for 
processing  complicated  signals  with  large  time-bandwidth  products,  because  of 
the  simplicity  of  the  synthesis  procedure  for  the  transversal  filler.  Although  the 
implementations  described  in  this  paper  as  well  as  other  existing  UTFs  provide 
heretofore  unavailable  signal  processing  capability,  there  are  important  immedi- 
ate applications  for  which  greater  capacity  is  required.  Microwave  acoustic  tech- 
nology, which  has  provided  the  basis  for  many  of  implementations  given  here, 
also  shows  great  promise  for  the  implementation  of  filters  to  meet  present  and 
future  advanced  requirements.  Because  of  the  close  relationship  between  the 
theories  and  techniques  employed  in  the  fields  of  acoustic  microwaves  and  electro- 
magnetic microwaves,  increased  collaboration  among  the  workers  in  these  two 
fields  will  make  significant  contributions  in  the  development  of  those  advanced 
devices. 
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